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The  performance  of  a modal  multidither  COAT  system  has  been  found 
to  be  quite  similar  to  an  equivalent  zonal  system  with  respect  to  target 
effects.  However,  the  whole-aperture  control  of  the  modal  system  results 
in  both  design  complications  and  alternatives  that  are  not  apparent  for 
zonal  systems.  These  are  especially  important  in  scenarios  where  low- 
order  aberration  correction  is  the  primary  requirement. 

All  multidither  systems  are  influenced  by  the  following  target 
effects: 

1.  When  the  target  spectrum  overlaps  the  dither  channels 

2 -1 

with  a normalized  power  density  Sc  > (G  Af)  , the 
affected  channels  will  be  dephased  (G  is  the  loop  gain 
and  Af  is  the  noise  equivalent  bandwidth) . 

2.  Large  targets  introduce  a dither  reduction  factor 
(DRF)  which  results  in  a non-optimum  gain  causing  de- 
layed convergence  and  increased  residual  uncompensated 
phase  error  at  convergence. 

A modal  transmitter  differs  from  a jonal  in  that: 

1.  Whole-aperture  control  does  not  dictate  the  loop  gain 
Increase  in  proportion  to  the  number  of  channels. 

2.  The  lower  order  modes  are  much  more  important  in  tur- 
bulence compensation  than  the  higher  order  modes,  whereas 
in  a zonal  system  essentially  all  channels  are  equiva- 
lent. 

3.  The  DRF  is  different  for  each  mode. 

4.  The  maximum  return  from  an  extended  target  may  occur 
at  a point  other  than  a converged  beam.  However,  the 
observed  effects  have  been  small  (i.e.,  ~ X/10  peak 
offset) . 
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5.  Dither  amplitudes  must  be  different  for  each  mode 
in  order  to  achieve  the  same  modulation  level. 

The  non-equivalence  of  each  channel  in  a modal  system  does  create  addi- 
tional design  complications.  However,  the  unequal  correcting  power  and 
whole-aperture  control  does  present  interesting  design  trade-offs  for 
adaptive  gain  and  channel  allocation  in  scenarios  with  target-limited 
performance. 

In  summary:  results  of  this  study  indicate  a modal  transmitter  will 

be  equivalent  to  a zonal  transmitter  in  the  presence  of  target  effects  (and 
perhaps  more  complex)  when  a large  number  of  channels  is  required  to  meet 
mission  requirements.  However,  in  scenarios  where  low-order  correction  is 
adequate  (i.e.,  tilt,  focus,  astigmatism)  a modal  system  can  be  conceived 
which  is  fundamentally  superior  and  less  complex  than  a zonal  equivalent. 
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INTRODUCTION 


1.1  BACKGROUND 

Technical  data  generated  by  the  COAT-Modal-Zonal  Comparison  Pregram, 
Contract  N60921-76-C-0122,  sponsored  by  the  Naval  Surface  Weapons  Center, 
monitored  by  Mr.  W.  P.  Altman,  NSWC,  is  presented  in  this  report.  The 
objective  of  this  program  was  to  assess  the  effects  of  target-COAT  inter- 
actions comparing  modal  multidither  with  zonal  multidither  and  phase-conju- 
gate systems.  This  program  paralleled  efforts  with  the  on-going  RADC/ARPA 
Contract  F30602-75-C-0086,^  which  is  supporting  zonal  multidither  and  phase 
conjugate  studies. 

A potential  problem  exists  when  target  modulation  overlaps  the  COAT 
sensing  bandwidth,  causing  erroneous  control  signals  to  be  generated.  Dif- 
ferences in  modal  and  zonal  transmitters  are  of  potential  importance  in 
closed-loop  performance  with  distributed  targets. 

1.2  DESCRIPTION  OF  COAT  OPERATION 

A COAT  system  uses  closed-loop  feedback  control  to  compensate  for 
time-varying  phase  distortions  which  would  otherwise  reduce  the  maximum 
achievable  energy  density  at  the  aimpoint.  These  systems  are  also  capable 
of  fine  pointing  and  tracking,  but  require  a strong  return  from  a highly 
localized  area  on  the  target's  surface  near  the  optical  axis  for  optimum 
performance.  A block  diagram  of  the  COAT  process  is  shown  in  Fig.  1-1. 

Two  different  sensing  techniques  are  considered  here,  one  based  on  energy 
detection  (outgoing  wave  multidither)  and  the  other  on  heterodyne  detection 
(return-wave  phase  conjugate). 

In  a multidither  system,  each  transmitter  element  (or  vibrant  struc- 
tural mode)  is  dithered  at  a distinct  frequency  to  phase  modulate  the  trans- 
mitted beam.  The  radiation  scattered  from  the  target  is  energy  detected; 
then  synchronous  detection  at  each  dither  frequency  is  used  to  recover  the 

*J.  Radley  et  al. , COAT-Target  Interactive  Modulation  Effects,  GRC  Report 
571-02-TR,  August  1976. 
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TRANSMITTER  ELEMENTS 


Figure  1-1,  Block  Diagram  of  COAT  Process 
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error  signal  corresponding  to  each  transmitter  element.  A "hill  climbing" 
algorithm  is  used  to  maximize  the  received  power  and  hence  the  power  on  the 
target. 

In  a phase  conjugate  system,  an  array  of  heterodyne  receivers,  which 
are  coincident  with  the  transmitter  elements,  measures  the  spatial  phase 
distribution  of  the  returned  wave  from  the  target  and  applies  the  negative 
of  this  to  the  transmitted  beam.  Assuming  reciprocity  of  the  phase  dis- 
turbances, the  transmitted  wavefront  cancels  them  and  produces  a diffrac- 
tion-limited beam  on  target. 

1.3  TARGET  EFFECTS 

Both  systems  have  been  shown,  both  experimentally  and  by  computer 
simulation,  to  perform  well  with  laser  wavefront  errors  and  atmospheric 
turbulence  for  point  targets.  However,  real  targets  are  distributed  and 
not  optically  smooth,  and  therefore  produce  spatially  complex  interference 
in  the  scattered  radiation  when  coherently  illuminated.  Target  rotation 
causes  this  spatially  varying  pattern  to  be  swept  across  the  receiving 
aperture (s)  in  time,  resulting  in  temporal  modulation  of  the  received  signal 
in  both  intensity  and  phase. 

In  a multidither  system,  target  effects  may  be  placed  into  two  cate- 
gories: temporal  and  spatial.  The  temporal  effects  depend  upon  the  size 
(D)  of  the  illuminated  area  on  the  target,  its  effective  rotation  rate  (ft) 
and  the  laser  wavelength  (A).  In  terms  of  power  spectrum,  the  target  modu- 
lation (TM)  noise  covers  frequencies  up  to  2ftD/A.  For  the  high  dither  fre- 
quencies the  TM  spectrum  may  not  overlap  the  sensing  bands,  but  the  analysis 
and  simulations  presented  here  show  that  there  is  noise  centered  about  each 
dither  frequency  which  reduces  SNR  and  increases  convergence  time.  If  the 
TM  spectrum  overlaps  all  the  dither  bands,  the  system  can  fail  to  converge 
altogether. 

Spatial  effects  depend  upon  the  brightness  distribution  of  the  target, 
which  in  turn  depends  on  the  target  geometry,  aspect  angle,  material 
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reflectance,  and  surface  roughness.  Adaptive  optics  (A/0)  systems  perform 
optimally  against  highly  localized,  single-glint  targets  because  these  give 
the  largest  dither  modulations  of  the  received  signal.  If  the  size  of  the 
glint  increases,  the  dither  signal  level  is  reduced  and  system  performance 
degrades.  This  effect  occurs  independent  of  the  target  dynamic  motion. 

In  a phase  conjugate  system,  an  extended  target  produces  a phase  in 
the  receiver  plane  which  has  a random  distribution  unrelated  to  the  desired 
correction  phase  and  a spatial  scale  the  same  as  that  of  the  intensity 
speckles.  Target  motion  produces  both  a translational  Doppler  frequency 
shift  from  the  optical  carrier  and  a frequency  spread  around  this  caused  by 
the  effective  rotation  of  the  target.  These  produce  temporal  variation  of 
the  phase  in  each  detector.  A target  which  has  a glint  of  large  spatial 
extent  causes  large  random  phase  errors  and  degrades  system  performance. 

1 . 4 STUDY  REQUIREMENTS 

The  Modal-Zonal  Comparison  Program  addressed  the  effect  of  targets  on 
COAT  systems  through  both  open-loop  and  closed-loop  simulation.  The  scope 
of  the  program  consisted  of  three  major  tasks: 

1.  Develop  the  formalism  to  describe  the  vibrational  modes  in  terms 
of  Zernike  polynomials,  which  would  be  compatible  with  computer 
implementation;  simulate  the  modal  transmitter  with  the  proper 
algorithms  to  permit  both  open-loop  and  closed-loop  operation; 
and  incorporate  the  MOWACS  (Modal  Wavefront  Control  System)^- 
modes  in  the  simulator. 

2.  Examine  the  open -loop  dither  error  signals  against  a point  target 
and  a variety  of  complex  targets,  to  determine  fundamental  differ- 
ences between  the  modal  and  zonal  transmitters. 

3.  Evaluate  illustrative  closed-loop  performance  against  a point 
target  and  a variety  of  complex  targets  with  and  without  a phase 
screen  atmosphere,  and  apply  parametric  analysis  to  performance 
predictions. 

Stodal  Wavefront  Control  Systems  (MOWArS)  Phase  II  Interim  Report , C.  Neufeld, 
Perkin-Elraer  Corporation  Report  12700,  September  1975. 
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1.5  SELECTION  OF  TARGETS  AND  SCENARIOS 


The  particular  choice  of  targets  and  scenarios  for  this  study  is 
important  because  it  serves  as  a basis  for  comparing  the  three  COAT  systems 


under  consideration.  It  is  apparent  from  the  variety  of  scenarios  en 
. , 1 . . . . . 


visioned  for  COAT  systems  that  undertaking  a parametric  study  of  perform- 
ance against  all  combinations  of  targets  and  scenarios  is  a monumental  (and 
unnecessary)  effort.  Much  of  the  interest  is  in  the  tactical  encounter  for 
which  the  atmospheric  effects  can  be  severe  and  hence  adaptive  compensation 
can  make  the  biggest  performance  improvement.  For  these  cases  atmospheric 
absorption  necessarily  limits  useful  ranges  to  less  than  20  km,  so  that 
usually  only  a portion  of  the  target  is  illuminated.  In  addition,  the 
target  aspect  angle  is  often  near  nose-on  (±30°)  for  many  tactical  scenarios, 
and  hence  can  be  represented  as  a simple,  single  dominant  feature,  as  shown 
in  Fig.  1-2.  In  this  way,  COAT  performance  can  be  better  quantified  in 
terms  of  the  effective  glint  size  on  the  target,  which  depends  on  the  sur- 
face radius  of  curvature  and  its  material  reflectance  properties. 


This  does  not  mean  that  all  targets  can  be  adequately  represented  by 
simple  sphere  segments  and  that  it  is  unnecessary  to  accurately  model  the 
particular  complex  target  of  interest.  However,  the  goal  of  this  study  is 
of  a more  fundamental  nature  and  is  meant  to  be  of  greater  applicability 
than  to  a particular  class  of  targets.  In  addition,  if  a COAT  system  does 
not  perform  well  on  a simple  sphere  segment  having  a single  maximum  radiance 
located  on  the  optical  axis,  then  it  almost  certainly  will  not  perform  well 
on  a complex  target  having  multiple  features. 


The  four  single-feature  targets  used  in  this  study  are  shown  in 
Fig.  1-3.  The  distributed  radiance  of  each  target  is  shown  on  a 14-level 
gray  scale.  The  display  scales  are  linear,  and  much  of  targets  A and  C are 
a factor  of  20  below  the  peak  value  and  hence  are  not  visible.  The  sizes 
and  radii  of  curvature  (0.2,  1,  and  4 m)  selected  roughly  correspond  to  a 
nose,  fuselage,  and  wing  section  of  an  aircraft  or  missile,  and  demonstrate 


J.  Radley  et  al.,  op.  cit. 


Distributed  Radiance  of  Single-Feature  Test  Targets  for  10.6-tim  Uniform  Illumination 
Display  is  14-level  linear  gray  scale.  Target  surface  is  aircraft  aluminum. 


the  dependence  of  effective  glint  size  on  target  geometry  and  reflectance 
(cf.  Fig.  1-3).  It  is  apparent  that,  of  the  four  targets,  target  A most 
closely  resembles  a point  target  against  which  the  COAT  system  should  per- 
form well.  Target  D is  the  flattest,  provides  the  least  spatial  discrim- 
inant for  a COAT  system,  and  should  cause  the  worst  performance. 

1.6  REPORT  ORGANIZATION 

The  remainder  of  this  report  contains  a general  description  of  target 
effects  on  COAT  systems  and  specific  comparisons  of  modal  and  zonal  multi- 
dither approaches.  Section  2 summarizes  the  nature  of  target  influence  on 
dither  system  performance.  Modal  and  zonal  system  characteristics  are 
compared  in  Sec.  3.  Section  4 provides  a description  of  the  simulation 
parameters  utilized  for  the  performance  evaluation  presented  in  Sec.  5. 

Additional  data  on  the  modal  simulation,  atmospheric  model,  and  aper- 
ture averaging  is  contained  in  the  Appendixes. 
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SUMMARY  OF  TARGET  EFFECTS  ON  DITHER  SYSTEMS 


The  spatial  extent  of  targets  and  the  temporal  variations  in  their 
backscatter  patterns  (l.e.,  speckle)  can  influence  COAT  performance.  These 
effects  have  been  studied  in  detail*  elsewhere.  The  results  are  summarized 
below  to  form  a basis  for  comparing  modal  and  zonal  COAT  characteristics. 

Real  dynamic  targets  must  be  characterized  by  both  spatial  and  tem- 
poral backscatter  properties.  These  properties  may  be  related  to  the 
spatial  and  temporal  resolving  power  of  the  COAT  system  to  determine  their 
Influence  on  convergence  properties. 

2.1  SPATIAL  EFFECTS  - DITHER  REDUCTION  FACTOR 

The  system  gain  and  dither  amplitudes  are  selected  in  order  to  achieve 
the  desired  convergence  time  and  the  final  level  of  convergence.  Two  first- 
order  design  considerations  are:  (1)  a hill-climbing  servo  can  approach 

the  maximum  to  within  a small  increment  (e) , and  as  the  gain  is  increased  £ 
is  reduced;  (2)  the  convergence  time  may  also  be  reduced  if  the  gain-band- 
width product  is  increased.  Thus  the  dither  level  is  selected  in  concert 
with  gain  parameters  in  order  to  achieve  the  maximum  gain  with  stable 
operation. 

The  finite  extent  of  a target  will  influence  the  depth  of  dicher  modu- 
lation sensed  by  the  COAT  system.  As  target  size  is  increased  irom  that  of 
a point,  the  strength  of  the  dither  is  modified  such  that  ;”.he  gain  is  no 
longer  optimal.  In  general  the  modulation  level  decreases  as  target  size 
Increases.  This  influences  the  system  performance  exactly  in  the  same  man- 
ner as  reduced  gain  which  delays  convergence  and  increases  the  residual  un- 
compensated phase  error  £.  In  most  instances  the  delay  rather  than  the 
residual  error  is  of  greatest  practical  concern. 

In  the  case  of  a zonal  transmitter  the  DRF  (Dither  Reduction  Factor) 
is  similar  for  each  transmitter  element.  It  becomes  of  significance  when 

*J.  Radley  et  al. , op.  cit. 


the  target  extent  approaches  that  of  the  single  element  diffraction  limited 
spot  size.  Thus,  as  the  number  of  channels  increases  (for  a fixed  total 
aperture),  the  DRF  is  less  important. 

The  effect  is  the  same  for  a modal  transmitter;  however,  its  magni-  ‘ 
tude  is  different  for  each  control  channel  since  the  diffraction  pattern 
of  each  mode  is  different.  In  general  the  lower  the  spatial  frequency 
content  of  the  mode  the  greater  the  DRF.  However,  as  is  shown  in  Sec.  3, 
the  relative  impact  on  the  nine  MOWACS  modes  (all  of  which  have  low  spatial 
frequency)  depends  on  the  specific  target  distribution. 

2.2  TEMPORAL  EFFECTS  - TARGET  MODULATION 

Temporal  effects  are  essentially  identical  for  modal  and  zonal 
systems. 

2.2.1  Target  Backscatter  Modulation 

The  complex  wavefront  emanating  from  a coherently  illuminated  target 
is  the  result  of  both  the  macroscopic  surface  geometry  which  is  much  larger 
than  an  optical  wavelength  and  the  microscopic  surface  which  may  be  on  the 
order  of  an  optical  wavelength  or  slightly  larger.  The  surface  microstruc- 
ture produces  random  phase  and  intensity  variation  in  the  backscatter  which 
is  noise  to  COAT  systems.  The  phase  and  intensity  profiles  in  the  receiver 
plane  are  illustrated  in  Fig.  2-1  for  two  different  target  sizes.  As  antic- 
ipated, the  larger  target  yields  higher  spatial  frequencies  in  the  back- 
scatter patterns. 

The  spatial  power  spectral  density  (PSD)  S (f  , f ) associated  with 

lx  y 

the  speckle  intensity  incident  on  the  receiving  aperture  is  related  to  the 
target  brightness  distribution  I(u,  v)  and  its  mean  <I>  by 

SI(fx*  V = <1>2t^fx*  V + ri(XRfx’  XRfy)]  (2_1) 

XJ.  W.  Goodman,  "Statistical  Properties  of  Laser  Speckle  Patterns,"  in 
Laser  Speckle  and  Related  Phenomena,  Topics  in  Applied  Physics,  Vol.  9, 

J.  C.  Dainty,  ed.  (Springer-Verlag,  New  York,  1975). 
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Figure  2-1.  Receiver  Plane  Intensity  (I/<I>)  and  Phase  (d>)  Corre' 
sponding  to  Two  Uniformly  Illuminated  Sphere  Segment 
Targets  Calculated  Using  Point  Scatterer  Model  (PSM) 
Computer  Program.  Range  ■ 5 km,  X ■ 10.6  ym;  phase 
scale  is  modulo  2ir.  (a)  Target  B:  diameter  * 10  cm 

radius  of  curvature  ■ 1 m;  (b)  Target  C:  diameter  * 

40  cm,  radius  of  curvature  ■ 1 m. 


where  TT(XRf  , XRf  ) is  the  autocorrelation  function  of  the  target  bright- 
i x y 

ness,  X is  the  wavelength,  and  R is  the  range.  The  spatial  frequencies 

f and  f correspond  to  the  x and  y directions,  respectively,  in  the  aper- 
x y 

ture  plane. 


The  detector  is  assumed  to  be  an  aperture-bounded  region  perpendic- 
ular to  the  line  of  sight,  and  produces  an  output  signal  proportional  to 
the  integrated  speckle  pattern  intensity  across  the  aperture  (i.e.,  a light 
bucket).  For  such  a detector  the  aperture  behaves  much  like  a two- 
dimensional  spatial  filter  which  attenuates  variations  at  different  spatial 
frequencies  to  different  degrees. 


The  spatial  PSD  Sn ( f , f ) of  the  speckle  intensity  that  is  passed 
u x y ^ 

by  the  aperture  is  related  to  ST(f  , f ) by  the  following  expression: 

J.  x y 


Sn(f  , f ) = ST(f  , f )|H(f  , f ) 
0 x y I x y 1 x y 


(2-2) 


where  H(f  , f ),  the  aperture  transfer  function,  is  the  two-dimensional 
x y 

spatial  Fourier  transform  of  the  pupil  function  of  the  receiver  aperture. 


Relative  motion  between  target  and  aperture  (e.g.,  a rotating  target) 
causes  the  speckle  pattern  to  be  swept  across  the  receiving  aperture  in  a 
direction  perpendicular  to  the  target's  effective  rotation  axis,  resulting 
in  a temporal-domain  PSD  Sg(f).  Since  all  spatial-frequency  components 
sharing  a common  projected  magnitude  along  the  direction  perpendicular  to 
the  effective  rotation  axis  give  rise  to  temporal  variations  at  the  same 
frequency,  the  PSD  of  the  detector  output  is  given  by  a one-dimensional 
integral  of  Sg(fx,  f^) . For  rotation  about  the  y-axis,  the  expression  for 
SQ(f)  is 

Srt(f)  - (v  r1  f Sn(f  , f ) df 
0 x 0 x y y 


f -f/v 


(2-3) 


Papoulis,  The  Fourier  Integral  and  Its  Applications,  (McGraw-Hill, 
New  York,  1962). 
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where  v - 2R ft  is  the  relative  velocity  in  the  x direction  at  the  aperture 
and  fl  is  the  target's  effective  rotation  rate  (i.e.,  the  projection  of  the. 
target's  rotation  vector  perpendicular  to  the  line  of  sight). 


2.2.2  Noise  Impact  on  COAT  Beam  Forming 

The  phasing  of  M transmitter  elements  results  in  an  on-axis  normalized 
intensity  in  the  target  plane  (Strehl  ratio)  given  by 


- (l/M2)  < f L 
- Ld“1 


23  exP 


m 


M 

23  exp 

Ln-l 


(2-4) 


where  is  the  target-plane  phase  of  the  field  from  the  ith  transmitter 

element.  If  we  assume  a Gaussian-distributed  phase  error  in  each  trans- 

2 

mitter  channel  of  variance  a,,  the  resulting  Strehl  ratio  is 


<I/1Q>  « exp  -02  + (1/M) (1  - exp  -cr2) 


(2-5) 


Since  it  can  be  shown  that  the  Strehl  ratio  is  not  highly  sensitive  to  the 
distribution  of  the  phase  errors,  the  Gaussian  assumption  detracts  little 
from  the  generality  of  the  following  results. 


Multidither  COAT  with  Noise 

For  a multidither  (M-D)  system,  such  as  that  shown  schematically  in 
Fig.  2-2,  the  variance  in  each  channel  is  produced  by  at  least  three  sources 
(1)  dither,  (2)  residual  phase  errors,  and  (3)  target  noise.  Thus  the  total 
variance  may  be  represented  as  follows: 

O2  - (B2/2)  + (a2/ G2)  + S0(fi)  • G2  • Af  (2-6) 

2 

where  0 is  the  dither  modulation  index,  ai  is  the  variance  of  the  initial 
phase  error  being  compensated,-  Af  is  the  noise-equivalent  bandwidth  (NEB) 
of  each  channel,  and  G^  and  G£  are  the  system  loop  and  electronic  gains  of 
the  control  channels  respectively. 
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Equations  2-5  and  2-6  may  be  combined  to  yield  the  steady-state  per- 
formance of  a multidither  system.  Table  2-1  indicates  the  base-line  par- 
ameters for  the  system.  The  resultant  performance  against  a simple  uniform 
square  test  target  is  illustrated  in  Fig.  2-3.  The  target  rotation  rate  is 
varied  to  scale  the  maximum  modulation  frequency  (f  ) in  the  spectrum. 

luoX 

For  this  example  the  target  is  assumed  to  be  very  small  such  that  aperture 
averaging  is  insignificant.  Distinct  regions  of  performance  which  approach 
nearly  perfect  phasing  and  total  random  phasing  are  evident.  The  transi- 
tion occurs  where  the  target  spectrum  SQ(f)  overlaps  the  sensing  frequencies 
with  high  power  spectral  density.  If  we  define  the  critical  level  of  noise 
in  a channel  as  that  amount  such  that  the  exponential  in  Eq.  2-5  drops  to 
e-1,  then  the  critical  noise  PSD  to  produce  channel  dropout  is 

S = 1/  (G2Af)  • (2-7) 

c 


Thus,  by  calculating  the  sensed  target  modulation  spectrum  it  is  possible 
to  estimate  the  dither  channels  which  will  be  dephased  by  the  target. 
Aperture  averaging,  gain,  and  bandwidth  trade-offs  may  be  employed  to 
minimize  target  influence  on  performance. 


Initial  Condition  Threshold 

The  previous  analytic  treatments  considered  the  targets  to  be  fully 
illuminated.  This  is  a good  approximation  for  small  targets  and  when  the 
beam  is  not  well  converged.  For  targets  which  are  physically  larger  than 
the  converged  spot  size,  the  effective  target  size  will  change  as  conver- 
gence is  approached.  This  results  in  a potential  threshold  process. 

With  the  beam  initially  unconverged  the  target  influence  may  be 
adequate  to  prevent  convergence.  However,  if  the  system  is  initially  con- 
verged the  effective  target  size  is  smaller  than  the  target's  physical  size. 
In  this  case  the  ability  to  maintain  convergence  may  be  dramatically  im- 
proved. With  intermediate  initial  conditions  the  system  will  either  decay 
to  a random  phase  condition  or  converge  depending  on  the  threshold. 


- 
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TABLE  2-1 

BASE-LINE  PARAMETERS 

8 *»  ±20°  = dither  modulation  index 

= n// 3 rad  » initial  RMS  phase  error 

Af  = 10ir,  30tt,  IOOtt  Hz  = noise  equivalent  bandwidth  (NEB) 
of  each  channel 

G^,  Gg  = system  loop  and  electronic  gains:  selected  to  main- 

tain a constant  gain-bandwidth  product  as  M and  Af 
are  varied. 

-3 

Tc  = 10  sec  =>  convergence  time  with  a point  target 

fp  = dither  frequencies:  sine-cosine  pairs  starting  at 

8 kHz  with  1.5  kHz  pair  spacing 

M = 10,  30,  100  = number  of  control  channels 

R = 5 km  * range  to  target 

A = 10.6  ym 

dT  = 60  cm  = transmitter  diameter 


Figure  2-3.  Strehl  Ratio  (I/Ig)  vs.  Maximum  Target  Modulation  Frequency 
(^max)  ^or  M-D  System  with  Uniform  Square  Target.  System 
noise  equivalent  bandwidth  (Af)  is  varied,  and  gain  is  ad- 
justed to  give  constant  convergence  time  (—1  msec). 

(a)  Af  - IOit  Hz;  (b)  Af  - 30tt  Hz;  (c)  Af  * IOOtt  Hz. 
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2.3  SUMMARY  OF  MULTIDITHER  COAT  PERFORMANCE 

The  extreme  sensitivity  to  system  gain  and  bandwidth  parameters  has 
been  developed  analytically.  Thus,  defining  general  regions  of  good  per- 

■ 

formance  depends  on  the  specific  system  parameters.  However,  it  is  pos- 

. . » ! 

sible  to  map  these  regions  adequately  to  summarize  and  illustrate  the 
trends  in  performance  as  target  and  system  parameters  are  varied. 

| | 

Figure  2-4  presents  a summary  of  multidither  (M-D)  performance 

i 

trends.  Three  performance  regions  (convergence,  threshold,  and  failure) 
are  shown  as  a function  of  glint  size  and  the  bandwidth  of  the  target 
spectrum. 

Large  targets  do  not  provide  an  adequate  discriminant  for  a M-D  hill 
climbing  algorithm  and  thus  result  in  failure  to  converge  independent  of 
the  target  modulation  bandwidth.  Small  targets  do  provide  a good  discrimi- 
nant and  the  system  will  converge  except  when  target  modulation  overlaps 
the  dither  frequencies.  As  shown  previously  the  extent  of  the  failure  re- 
gion due  to  noise  overlapping  the  dithers  depends  upon  the  system  param- 
eters (cf.  Fig.  2-3). 

U m 

In  the  threshold  region,  M-D  performance  depends  upon  the  initial 
conditions:  level  of  convergence,  aimpoint,  etc.  If  the  system  is  at  or 

near  convergence  it  should  be  able  to  maintain  convergence  in  this  region. 

However,  if  the  initial  conditions  are  poor,  the  system  is  likely  to 
transistion  to  a failure  mode.  As  shown  in  Fig.  2-4,  this  region  is 
bounded  by  target  sizes  corresponding  to  the  converged  spot  size  and  single 
element  spot  size.  Within  these  bounds  the  effective  target  size  is  deter- 
mined by  the  illumination  and  thus  the  threshold  effect. 


COMPARISON  OF  SYSTEM  CHARACTERISTICS 

Both  modal  and  2onal  systems  considered  here  are  based  on  the  same 
signal  processing  and  convergence  algorithm.  In  this  section  the  differ- 
ences in  transmitter  characteristics  and  any  restrictions  they  imposed 
on  the  parameter  selection  for  the  signal  processing  are  briefly  noted. 

3.1  ATMOSPHERIC  COMPENSATION 

Other  studies  have  shown  quantitatively  that  for  low-order  correc- 
tions (tilt,  focus,  etc.)  a modal  correction  can  be  tailored  to  provide 
compensation  with  fewer  channels  than  a zonal  system.  However,  when  higher 
order  aberrations  must  be  compensated  the  modal  advantage  is  lost  entirely. 

3.2  BEAM  CONTROL/DITHER 

The  differences  in  modal  and  zonal  concepts  are  clearly  apparent  in 
the  beam  control  characteristics.  For  an  M channel  system  each  zonal  ele- 
ment controls  1/M  of  the  total  power  with  an  influence  function  in  the 
target  plane  corresponding  to  the  single-element  diffraction  pattern.  In 
the  modal  case  each  mode  perturbs  the  entire  transmitted  beam  with  the 
diffraction  pattern  corresponding  to  the  mode  profile  as  the  influence 
function  in  the  target  plane.  Thus  each  zonal  channel  is  essentially 
identical  in  dither  and  beam  control  characteristics,  whereas  each  modal 
channel  is  different.  A modal  transmitter  will  have  the  additional 
complexity  of  (1)  different  excitation  amplitudes  for  each  mode  to  achieve 
desired  modulation  depth,  (2)  different  response  to  an  extended  target  in 
each  channel. 

3.3  BEAM  CONTROL  INFLUENCE  ON  ELECTRICAL  LOOP  DESIGN 

In  a zonal  transmitter  as  the  number  of  channels  is  increased  the 
loop  gain  must  be  increased  to  account  for  the  lower  signal.  In  the  modal 
case,  whole  aperture  control  in  principle  would  permit  the  gain  to  remain 
constant  as  more  channels  are  added.  However,  in  practice  the  residual 
modulation  at  convergence  (due  to  dither)  results  in  the  requirement  for 
reduced  dither  levels  as  channels  are  added  and  higher  gain,  just  as  in  the 
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4 SIMULATION  BASE-LINE  PARAMETERS 

As  shown  in  Fig.  4-1,  the  computer  simulator  used  in  this  study  con-  * 
sists  of  three  basic  functional  components:  (1)  target  simulator  (TARGSIM) 

which  calculates  the  target's  laser  radar  cross  section  based  on  its  geo- 
metrical description,  materials  reflectance  properties  and  motion  parameters; 
(2)  point  scatterer  model  (PSM)  which  calculates  the  backscattered  radiation 
from  the  target  propagated  to  the  receiver  plane;  and  (3)  COAT  simulator, 
which  includes  the  receiver,  signal  processor,  and  adaptive  transmitter. 

The  following  COAT  simulators  are  used  for  the  closed-loop  comparisons 
presented  in  this  section:  MUDuATS  (MUltiDither  £OAT  Analysis  ^Through 

Simulation),  TOPCATS  (The  Optical  £hase  Conjugate  Adaptive  Techniques  Simu- 
lation) , and  DEMCATS  (DEformable  Mirror  £OAT  Analysis  through  Simulation) • 

The  first  two  are  described  in  complete  detail  by  Radley  et  al.  and  are  only 
shown  schematically  here  (in  Figs.  4-2  and  4-3)  to  indicate  how  each  is  con- 
structed. Some  of  the  zonal  multidither  and  phase  conjugate  closed-loop 
simulations  presented  in  the  reference  are  also  presented  in  Sec.  5 for  com- 
parison with  the  modal  multidither  simulations  accomplished  for  this  study. 

The  modal  simulator  (DEMCATS)  is  described  in  detail  in  Sec.  4.1,  which 
also  includes  the  rationale  for  selecting  the  system  parameters. 

4.1  MODAL  SYSTEM  BASE-LINE 

The  computer  program  (DEMCATS),  which  was  constructed  to  model  a modal 
multidither  COAT  system,  is  shown  schematically  in  Fig.  4-4.  The  only  differ- 
ence between  it  and  the  zonal  simulator  (MUDCATS) , shown  in  Fig.  4-2,  is  in 
the  error  signal  processor  which  accomplishes  the  transformation  from  signal 
level  (dither  plus  error  signal)  per  channel  to  total  mirror  displacement 
(as  described  in  Appendix  A).  The  zonal  transmitter  can  be  represented  by 
a similar  matrix.  The  simulator  description  and  selection  of  parameters 
are  presented  in  the  following.  Additional  data  on  the  modal  base  line 
development  and  evaluation  are  provided  in  Appendix  B. 


Radley  et  al.,  op. 


cit. 
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PROGRAM  FLOW 


Scattering  Model/COAT  Simulator  Interface.  Receiver  is  depicted  schematically  only 
extended  and  multiple  receivers  are  included.  Also,  point  scatterers  are  not  neces- 
sarily in  a regular  array. 


r i 

Figure  4-3.  Block  Diagram  of  TOPCATS  (The  Optical  Phase  Conjugate 

Adaptive  Techniques  Simulation).  Ideal  optical  hetero- 
dyne receivers  and  phase  comparators  are  assumed. 


Figure  4-4.  Block  Diagram  of  DEMCATS  (DEfonnable  Mirror  COAT  Analysis  Through  Simulation) 


4.1.1  Description  of  Simulator 

The  computer  program  DEMCATS  has  the  following  features  and  capabil-  * 
ities.  It  is  fully  closed  loop  with  a complex  target  as  an  integral  part 
of  the  system  and  incorporates  all  of  the  important  signal  processing 
features:  energy  detection  with  various  aperture  sizes  and  shapes  (e.g., 

point,  annulus,  coincident,  mesh),  automatic  gain  control  with  signal  clip- 
ping, synchronous  detection,  and  bandpass  filtering.  The  number  of  modes 
(each  represented  as  a sum  of  Zernike  polynomials)  is  selectable,  as  are 
the  dither  frequencies,  interdither  spacings,  system  cycle  time,  dither 
levels,  system  gains,  and  filter  cutoff  frequencies.  Propagation  is 
accomplished  by  a direct  ray  trace. 

The  system  operates  in  the  following  manner.  Once  the  scenario 
parameters  and  point-scatterer  representation  of  the  target  have  been  de- 
termined, the  initial  phase  error  is  applied  to  the  system,  usually  by 
randomizing  the  phase  of  the  individual  transmitter  modes  from  -tt/2  to  tt/2, 
but  in  some  cases  by  generating  a phase  screen  atmosphere  with  the  atmos- 
pheric model  described  in  Appendix  C.  Perfect  focus  and  no  tilt  error  are 
assumed  so  that  each  mode  has  only  to  correct  its  phase  error  (except  in 
the  case  of  atmospheric  phase  error) . Each  mode  undergoes  both  a high- 
frequency  dither  (8  kHz  and  above)  and  a lower  frequency  correction  (-100- 
300  Hz). 

This  beam  is  propagated  to  the  target  and  back  to  the  receiver  (usually 
an  annulus  with  1.2:1  ratio  of  O.D.  to  I.D.)  which  is  assumed  to  have  zero 
noise.  This  signal  is  then  multiplied  by  the  detector  gain  (GD  in  Fig.  4-4) 
to  convert  target  radiance  to  volts.  The  automatic  gain  control  (AGO)  is  an 
analog  divider  which  low-pass  filters  a portion  of  the  received  signal  and 
divides  this  into  the  unfiltered  component.  This  is  necessary  for  stability 
considerations  for  large  numbers  of  channels  and  for  improved  performance 
in  the  presence  of  target  noise.  The  analysis  of  the  AGC  is  reported  else- 

i 

wh«*re.  * 


J.  Radley  et  al, , op. 
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I The  high-  and  low-pass  filters  provide  large  attenuation  at  fre- 

quencies outside  the  dither  frjquency  bands.  The  saturable  device  is  a 
"soft"  clipper  which  limits  the  signal  for  approximately  the  first  0.5  msec 
until  the  low  pass  filter  in  the  AGC  has  had  sufficient  time  to  charge  and 
provide  a nonzero  voltage.  Otherwise,  the  initial  instability  drives  one 
or  more  channels  to  a state  which  will  produce  only  a local  maximum  at 
"convergence."  The  clipping  mechanism  is  logarithmic  rather  than  a sharp 
cutoff  to  reduce  the  characteristic  high-frequency  ringing  which  would  ex- 
tend into  the  dither  bands. 

After  amplification  by  the  electronic  gain  (G^,) , the  signal  is  then 

synchronously  detected  at  each  of  the  dither  frequencies.  This  detection 

process  selects  only  the  signal  at  each  particular  dither  frequency  and 

shifts  this  to  zero  frequency.  The  error  signal  is  then  filtered  by  a 

series  of  five  low  pass  filters.  The  most  important  filter  is  LPF^  which 

has  a 10  Hz  break  point  and  eliminates  most  of  the  noise  which  is  at  higher 

frequencies  than  the  error  signal  (which  is  DC).  The  other  filters  have 

a substantially  higher  cutoff  frequency  and  reduce  intermodulation  noise. 

1 2 

The  particular  choice  of  these  is  discussed  in  detail  in  other  reports.  ’ 

The  dither  signal,  producing  a modulation  index  8,  is  added  to  the 
error  signal  applied  to  each  transmitter  element  in  order  to  tag  the  error 
signal  per  channel  prior  to  the  next  cycle.  The  convergence  algorithm  is 
a "hill  climbing"  type  which  seeks  to  minimize  the  error  signal  thereby 
maximizing  the  total  received  power. 

This  process  is  repeated  every  10  ysec  (the  simulated  cycle  time)  or 
less,  depending  upon  the  required  sampling  rate  for  dither  or  target  modu- 
lation frequencies. 


^W.  B.  Bridges  et  al..  Coherent  Optical  Adaptive  Techniques  (COAT) , Hughes 
Research  Laboratories,  Technical  Reports  No.  1 & 2,  July  and  October  1973. 

2 

J.  E.  Pearson  et  al..  Coherent  Optical  Adaptive  Techniques  (COAT) , Hughes 
Research  Laboratories,  Final  Technical  Report,  RADC-TR-75-46,  January  1975. 
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The  selection  of  system  parameters  for  all  closed  loop  simulations 
Is  based  on  the  following  criteria: 

a.  The  parameters  are  the  most  favorable  to  performance 
against  an  ideal  (point)  target,  for  which  a convergence 
level  of  at  least  90Z  must  be  attained  within  2 msec. 

b.  They  are  physically  realizable  using  present  or  antici- 
pated state-of-the-art  hardware. 


In  order  to  isolate  only  the  effects  of  target  noise  and  provide 
errors  which  the  systems  should  be  capable  of  correcting,  the  following 
conditions  are  incorporated  in  the  simulators:  infinite  detector  signal- 

to-noise  ratio,  no  tilt  error,  and  a lens  at  the  transmitter  focused  on  the 
target.  The  base-line  system  utilizes  the  first  9 MOWACS  modes.  The  over- 
all transmitter  diameter,  60  cm,  is  within  the  50-150  cm  range  proposed  for 
most  scenarios. 

The  base-line  initial  condition  is  random  peak  mode  amplitude  dis- 
tributed — rr/2  to  ir/2  in  phase  with  no  atmosphere.  This  represents  a phase 
error  which  each  mode  is  able  to  correct.  Because  the  RMS  phase  error  char- 
acterizes the  initial  level  of  convergence  and  because  only  nine  samples  of 
a uniform  distribution  can  give  significantly  different  RMS  on  each  run,  the 
same  initial  phases  were  used  for  the  base-line  system  simulations.  This 
phase  distribution  has  an  RMS  error  of  1.7  rad  and  a peak-to-peak  error  of 
6.28  rad. 


Refer  to  Fig.  4-4;  the  following  parameter  values  were  chosen: 
Number  of  channels  **  9 

Dither  frequencies  ■ 8-20  kHz,  Af  * 1.5  kHz 

Gains:  8 ■ 0.503  (±28.8°  dither  modulation  for  Mode  1) 

Gd  - 12,000 

G£  - 500 
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Filter  break  frequencies:  LPF^  * 10  Hz 

LPF2_5  = 5 kHz 
AGC  LPF  = 1 kHz 
2 HPF’s  = 1 kHz 
2 LPF's  = 50  kHz 

Saturable  device  cutoff  (logarithmic)  * 50  V 

4.1.2  Modulation  Index 

The  first  step  was  to  determine  the  dither  level.  In  a zonal  multi- 
dither system  this  is  ±20°  based  on  a trade-off  between  required  signal  level 
and  linearity  of  operation.  The  amount  of  control  of  on-axis  power  on  target 
is  shown  in  Sec.  5 for  7,  19,  and  100  channels  each  with  the  same  overall 
transmitter  size.  As  expected,  the  percent  modulation  decreases  with  the 
number  of  elements.  For  7 elements,  a ±20°  dither  at  90%  convergence  corre- 
sponds to  ±6.4%  power  modulation.  This  value  and  an  alternative  dither  level 
equivalent  to  the  19-channel  zonal  transmitter  was  sought  to  match  the  MOWACS 
transmitter.  To  compensate  for  a decrease  in  dither  level  for  the  19  channel 
case,  a proportional  increase  in  signal  gain  was  made. 

Each  MOWACS  mode  has  a different  spatial  structure  and  hence  produces 
a different  amount  of  target  plane  power  modulation  for  a given  amount  of 
dither.  To  make  a modal-zonal  comparison,  we  chose  the  amount  of  dither  for 
each  mode  to  give  the  same  signal  modulation  at  the  same  convergence  level. 
Each  mode  was  assigned  a different  peak  dither  amplitude  (±29°  for  mode  //l, 
±9°  for  mode  #4,  etc.)  so  each  would  possess  the  same  "optical"  gain,  i.e., 
the  same  power  modulation  on  a point  target.  The  dither  levels  in  mode 
amplitude  are  listed  in  Table  4-1. 

4.1.3  Dither  Frequencies 

Initially,  the  intent  was  to  use  the  actual  MOWACS  dither  frequencies 
(2.28-13.3  kHz);  however,  this  was  not  done  for  initial  comparisons  for  the 
following  reason.  More  target  noise  occurs  at  lower  frequencies  for  all 
target  rotation  rates,  hence  the  low  dither  channels  encounter  more  target 
noise.  Thus  the  same  frequencies  and  interdither  separations  as  the  zonal 
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TABLE  4-1 

DITHER  LEVELS  FOR  MODAL  TRANSMITTER 


Mode  No. 
1 
2 

3 

4 

5 

6 

7 

8 
9 


Dither  (tdeg) 
28.8 

13.5 

13.5 
9.2 
10.7 
12.1 
12.1 
15.3 

13.5 
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case  were  used  (f^  * 8-20  kHz,  Af  * 1.5  kHz).  This  permits  a direct  com- 
parison of  full  aperture  modal  beam  control  vs.  segmented  control.  The 
8-kHz  lowest  dither  frequency  creates  a region  where  low-frequency  target 
noise  does  not  overlap  the  dither  frequencies.  This  permits  a clearer 
illustration  of  the  target  effects  as  dynamic  rates  are  varied. 

4.1.4  Initial  Transmitter  Phase  Error 

The  next  step  was  to  determine  an  initial  phase  error  comparable  to 
that  assumed  in  the  zonal  case  (random  0 to  2u  per  element) . Also  the 
spatial  frequency  could  not  exceed  the  correction  capability  of  the  first 
nine  MOWACS  modes.  The  peak  mode  error  was  random  -tt/2  to  tt/2  and  a number 
of  randomizations  were  made  to  select  the  same  initial  RMS  error  as  the 
zonal  case  with  a nonzero  initial  error  for  each  mode.  However,  a major 
difference  between  modal  and  zonal  is  the  modulo  2ir  factor.  In  the  zonal 
case,  an  element  can  be  2tt  out  of  phase  with  the  other  elements  and  still 
produce  near-diffraction-limited  performance  because  each  is  independent 
(spatially)  of  the  others.  In  the  modal  case,  the  mode  errors  are  additive, 
which  can  lead  to  peak-to-peak  errors  too  large  for  the  system  to  correct. 
One  mode  2ir  out  of  phase  produces  a large  reduction  in  the  power  on  target. 
(See  Appendix  B,  Sec.  B.2,  for  sensitivity  to  phase  error.)  Figure  4-5 
illustrates  the  target  plane  beam  profiles  which  were  used  in  the  comparison 
between  modal,  zonal,  and  phase  conjugate  closed-loop  performance.  Final 

modal  runs  used  the  transmitter  with  <p  =6.28  rad. 

P-P 

4.1.5  Automatic  Gain  Control 

Zonal  results  also  demonstrated  the  importance  of  an  AGC  for  minimiz- 
ing system  gain  variation  resulting  from  the  low-frequency  signal  modulation 
caused  by  low  target  rotation  rates.  An  AGC  was  also  found  to  improve  zonal 
multidither  COAT  performance  against  targets  producing  higher  frequency 
noise.  Hence,  it  was  included  in  DEMCATS  prior  to  making  any  complex  target 
runs.  Basically,  the  AGC  divides  the  received  signal  by  a 1-kHz  low-pass- 
filtered  version  of  the  same  signal.  In  this  way,  signal  fluctuations  are 
normalized  to  the  mean  of  the  received  signal. 
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To  prevent  initial  instability  caused  by  large  gains  from  an  un- 
charged LPF  in  the  AGC,  the  error  signals  are  disabled  for  the  first  0.125 
msec.  Any  further  instabilities  incurred  during  the  engagement  due  to 
signal  dropout  (speckle  nulls)  are  checked  by  a "soft"  voltage  clipper. 

The  soft  clipper  is  logarithmic  for  voltages  above  the  clipping  level.  Its 
main  purpose  is  to  limit  large  gain  fluctuations  created  by  the  AGC  (during 
signal  dropout) . 

4.1.6  Aperture  Averaging 

Multidither  COAT  systems  (both  modal  and  zonal)  sense  the  intensity 
or  power  of  the  radiation  scattered  from  the  target  and  use  synchronous 
detection  to  selectively  measure  the  signal  modulation  at  each  of  the 
separate  beam  tagging  (or  dither)  frequencies.  Any  spurious  intensity 
modulation  at  these  frequencies  is  noise  which  competes  with  the  signals 
in  the  affected  channels. 

Most  real  targets  are  not  optically  smooth  and  hence  produce  a 
spatially  random  intensity  pattern  in  the  receiver  plane  which  can  be 
characterized  by  a probability  density  function  and  a power  spectral  den- 
sity (or  its  Fourier  transform,  the  autocorrelation  function).^-  Any  method 
which  reduces  the  dynamic  range  of  these  fluctuations  and  its  extent  in 
the  frequency  domain  reduces  its  impact  on  multidither  system  performance. 

Aperture  averaging  is  one  such  method  of  reducing  target  noise.  In 
fact,  a larger  receiver  aperture  is  also  beneficial  in  terms  of  increasing 
the  received  signal  level  with  respect  to  the  detector  and  system  electronic 
noise.  One  of  the  easiest  to  implement  from  an  optical  design  standpoint 
is  an  annular  receiver  which  is  coaxial  with  the  transmitter.  In  this  way, 
the  same  beam-expanding  telescope  can  be  used  to  transmit  and  receive.  A 
configuration  which  is  much  more  difficult  to  implement  is  a receiver  which 
is  coincident  with  the  transmitting  aperture.  In  this  case,  a high-power 
beam  splitter  is  required  and  the  optical  system  is  faced  with  some  of  the 

*J.  Radley  et  al.,  op.  cit. 
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same  problems  as  return-wave  phase  conjugate.  It  has  been  shown  elsewhere 
that  both  configurations  are  nearly  equal  in  terms  of  reducing  target  modu- 
lation noise. 

Rather  than  make  a large  number  of  closed-loop  simulations,  it  is 
more  instructive  and  considerably  more  efficient  to  use  the  techniques  for 
calculating  open-loop  aperture-averaged  PSD,  probability  density,  and  con- 
trast ratio  developed  in  Appendix  D to  determine  the  effect  of  aperture 
averaging  and  provide  an  initial  estimate  of  which  targets  and  rotation 
rates  are  potentially  most  stressful  to  system  performance. 

Briefly  summarized,  the  target  spatial  PSD  Is  the  sum  of  a DC  com- 
ponent and  higher  frequency  component  calculated  from  the  autocorrelation 
of  the  target  brightness  distribution.  It  requires  an  additional  step  if 
the  spatial  structure  of  the  adaptive  illumination  is  to  be  included,  but 
for  this  analysis  it  has  little  effect  and  uniform  illumination  is  assumed. 
(This  is  valid  provided  that  the  beam  is  larger  than  the  target,  which  is 
the  case  here  when  the  beam  is  unconverged . ) 

The  spatial  PSD  Is  then  multiplied  by  the  magnitude  squared  of  the 
aperture  transfer  function,  and  this  is  integrated  in  the  spatial  frequency 
component  which  corresponds  to  the  direction  parallel  to  the  target's  rota- 
tion axis.  Integrating  the  spatial  PSD  in  both  the  x and  y spatial  fre- 
quencies gives  the  number  of  correlation  cells  in  the  aperture  and  hence 
the  speckle  probability  distribution  and  contrast  ratio  (cf.  Appendix  D) . 

Finally,  it  is  scaled  by  2Rft  to  transform  it  into  the  temporal  fre- 
quency domain.  This  process  has  been  programmed  on  a digital  computer  and 
makes  extensive  use  of  fast  Fourier  transforms  so  that  parametric  variation 
of  target  geometry,  receiver  geometry,  and  scenario  parameters  (range  and 
wavelength)  is  fast  and  efficient.  Finally,  this  analytic  method  produces 
ensemble  average  results  without  the  necessity  of  generating  a large  number 
of  statistical  samples  to  adequately  predict  expected  performance. 

\j.  Radley  et  al.,  op.  cit. 
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The  following  receiver  geometry  is  considered  here:  70-cm  o.d., 

60-cm  i.d.  annulus.  Both  a one-dimensional  cut  and  the  y-spatial-f requency- 


integrated,  aperture-averaged  PSD's  are  shown  in  Fig.  4-6a.  Note  that 
the  latter  is  the  aperture-averaged  PSD  (in  the  spatial  frequency  domain) 
for  an  infinite,  uniform  target  rotating  about  the  y axis. 


1 1 

Li 


The  real  test  is  how  these  perform  in  terms  of  reducing  target  noise. 
Rather  than  using  uniform  square  targets  of  various  sizes  (whose  PSD's  are 
simply  triangular)  one  of  the  targets  used  for  the  open-loop  simulations  in 
Sec.  5 is  used  here.  This  target  is  an  intermediate  case  of  those  considered 
in  this  study.  The  aperture-averaged  PSD  is  shown  in  Fig.  4-6b. 


However,  the  PSD  is  only  half  the  story.  The  probability  density 
function  shows  which  values  of  intensity  are  most  probable,  and  the  ratio 
of  the  RMS  to  the  mean  is  defined  as  the  contrast  ratio.  Aperture  averaging 
may  produce  a significant  increase  in  the  mean  signal  level  and  reduction  in 
the  fluctuations.  The  probability  density  functions  are  shown  in  Fig.  4-7. 


In  order  to  scale  the  target  modulation  spectrum  (PSD)  presented  in 
Fig.  4-6b  to  the  temporal  frequency  domain  for  comparison  to  the  dither 
sensing  frequencies,  the  spatial  frequency  scale  must  be  multiplied  by  2Rf2 
and  the  vertical  scale  (PSD)  divided  by  2Rf)  (to  keep  the  total  integrated 
power  constant)  as  described  in  Appendix  D.  These  are  presented  in  Figs. 

4-8  and  4-9,  for  the  four  spherical  targets  considered  in  Sec.  5,  at  two 
different  rotation  rates.  It  is  apparent  from  these  that  there  is  more  tar- 
get noise  at  lower  frequencies  for  all  cases;  hence  operation  at  higher 
dither  frequencies  is  very  desirable  for  reducing  target  effects.  Also, 
when  the  adaptive  beam  is  initially  unconverged  and  illuminates  a larger 
area  of  the  target,  as  the  extended  target  cases  presented  here,  target  noise 
may  extend  to  tens  of  kilohertz  for  tenths  of  radians  per  second  rotation 
rates.  Hence,  a system  which  can  operate  with  fewer  control  channels  and  at 
higher  dither  frequencies  will  be  able  to  operate  above  the  target  modula- 
tion frequencies  encountered  in  many  engagement  scenarios.  A 60-100-channel 
zonal  multidither  system,  which  is  constrained  to  dither  frequencies  down 
to  a few  kilohertz,  will  suffer  more  severe  performance  degradation. 
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2-D  RCVR, 1-D  RCVR 


(a)  One-Dimensional  and  Two-Dimensional  PSDs  of  an  Annular  Receiver.  The 
latter  also  corresponds  to  the  target  modulation  spectrum  for  the 
idealized  case  of  a uniform,  infinite  target  rotating  about  the  y axis. 


I 

( 


( 
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<P$D>  (dB) 


f (m“^) 
x 

Original  Target  PSD,  Aperture  Averaged  PSD 

(Point  Receiver) 

(b)  Target  Modulation  Spectra  for  Two  Receiver  Geometries.  Parameters: 
sphere  segment  target  (dQ  * 40  cm,  rc  ■ 1 m)  uniformly  illuminated, 
range  ■ S km,  X ■ 10.6  ym. 


Figure  4-6.  Ensemble  Average  Power  Spectral  Densities  (PSDs)  Corre- 
sponding to  an  Annular  Receiver  (O.D.  * 70  cm,  I.D.  ■ 

60  cm). 


. 
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(a)  Speckle  probability  distribution  vs.  number  of  correlation 
cells  (M)  in  aperture. 


Aperture  Shape 

Diameter  (cm) 

M 

CR  = al/<I> 

Circular 

60 

7.58 

0.364 

Annular 

o.d.  = 85,  i.d.  = 60 

12.05 

0.288 

Annular 

o.d.  =70,  i.d.  =60 

9.47 

0.325 

Annular 

o.d.  = 61 , i.d.  = 60 

8.56 

0.342 

1 Computer  calculated  parameters  (M, 

CR)  for  each  aperture 

Figure  4-7.  Effect  of  Receiver  Geometry  On  Probability  Density  Function 
of  Direct  Detected  Signal.  Parameters  for  (b)  given  in 
Fig.  4-6b. 


D1 ther 
Frequencies 


Figure  4-8.  Ensemble  Average  Power  Spectra  of  Target  Modulation  (Noise)  for  Uniform  Illumination. 
Continuous  and  dashed  lines  Indicate  point  and  annular  receivers,  respectively. 

Target  rotation  rate  is  0.25  rad/sec. 


(9P)  <0$d>  L 0 * 0 (9P)  <QSd>  Lo*0 


Computer  run  time  requirements  constrained  us  to  use  only  the  unre- 
solved (dg  * 10  cm)  target  for  the  illustrative  closed  loop  simulations  in  . 
Sec.  5.  However,  it  is  clear  from  the  open  loop  PSD's  in  Figs.  4-8  and  4-9 
that  the  cases  where  more  surface  area  is  illuminated  or  where  the  glint  is 
less  localized  will  be  more  stressful  to  system  performance. 

4.1.7  Atmospheric  Turbulence 

In  Sec.  4.1.4  an  initial  set  of  random  amplitude  modal  perturbations 
is  applied  to  the  transmitter  surface  to  approximate  an  atmospheric  phase 
screen.  The  implication  is  that  such  a two-dimensional  perturbation  in 
the  transceiver  plane  accounts  for  most  of  the  phase  variations  caused  by 
propagation  through  the  turbulent  atmosphere.  Using  an  atmospheric  model 
described  in  Appendix  C,  two  different  atmospheric  phase  screens  were  gen- 
erated for  use  in  the  modal,  zonal,  and  phase  conjugate  COAT  simulators. 

The  model  parameters  were  assigned  the  following  values: 

1.  Outer  scale,  Lq  * 10  m 

2.  Effective  width  of  the  phase  perturbing  region,  d = 100  m 

2 -13  -2/3 

3.  C m 10  (m  ) - moderate  turbulence 

n 

4.  * 10  ^ (m  - weak  turbulence 

n 

Using  the  Fig.  C-l,  or  Eq.  C-l,  the  RMS  phase  fluctuations  (a^  values)  for 
the  moderate  and  weak  turbulence  are  1.19  rad  and  0.376  rad,  respectively. 

The  modal  transmitter  point  locations  were  assigned  phase  screen  values 
that  corresponded  to  their  physical  location.  As  the  modal  correction  is 
engaged,  the  phase  screen  is  compensated  by  the  negative  wavefront  of  the 
transmitter.  Using  the  same  atmospheric  phase  screen  for  the  zonal  trans- 
mitter (and  phase  conjugate)  requires  sampling  the  phase  screen  at  each  of 
the  transmitter  element  centers  and  determining  the  x and  y phase  slope 
across  each  of  the  elements.  A simple  linear  approximation  using  two  phase 
screen  points  for  each  orthogonal  direction  about  each  transmitter  element 
center  is  employed.  The  mean  phase  slopes  are  applied  to  the  diffraction 


pattern  of  each  transmitter  element,  but  cannot  be  directly  compensated 
for  because  of  the  single-direction,  piston-like  movement  of  the  trans- 
mitter elements. 

Figure  4-10  illustrates  the  atmospheric  phase  screen  used  in  the  COAT 
simulators.  The  values  are  normalized  to  the  maximum  phase  excursion.  The 
minus  sign  over  some  of  the  array  values  indicates  negative  phase  values. 
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Figure  4-10.  Normalized  Atmospheric  Phase 
Screen.  Parameters:  Ln  ■ 10 


Screen.  Parameters:  Lq  ■ 10  m* 

d - 100  m.  For  C2  - 10-14  m_2/3, 
maximum  phase  excursion  <■  0.832 
rad;  for  c£  - 10~13  m-2/3,  $mair  . 
2.631  rad.  Minus  signs  indicate 
negative  phase  values. 
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5 PERFORMANCE  EVALUATION 

Utilizing  the  simulation  techniques  described  in  Sec.  4 both  open- 
and  closed-loop  evaluations  are  considered  in  this  section  to  investigate 
the  factors  identified  in  Sec.  3.  In  addition  to  the  modal  and  zonal 
multidither  simulations,  phase  conjugate  simulations  are  also  provided  to 
complete  the  COAT  comparisons. 


5.1  OPEN-LOOP  M-D  COMPARISONS 

There  are  several  important  things  to  note  from  the  open  loop  evalu- 
ations which  follow:  (1)  each  mode  has  a different  optical  gain  (i.e., 

for  the  same  peak  dither  amplitude  with  a point  target  a different  modula- 
tion index  results);  (2)  the  ranking  of  modes  by  optical  gain  depends  on 
the  particular  radiance  distribution  of  the  target;  (3)  for  all  of  the 
distributed  targets,  some  of  the  modes  produce  a maximum  in  the  integrated 
radiance  (and  hence  received  signal)  at  nonzero  mode  amplitude;  and  (4) 
these  three  effects  are  different  than  for  a zonal  transmitter. 

5.1.1  RMS  Phase  Error  of  Transmitter  Surface 

The  fundamental  difference  between  the  modal  and  zonal  transmitter 
is  the  control  of  the  transmitter  surface.  Each  control  channel  for  the 
zonal  transmitter  affects  a segment  of  the  total  laser  energy  incident 
across  its  aperture,  whereas,  in  the  case  of  the  modal  transmitter,  the 
entire  surface  is  controlled  by  each  of  the  control  channels.  To  deter- 
mine the  relative  control  over  the  transmitter  surface,  we  can  compare 
transmitter  surface  RMS  vs.  peak  phase  excursion  plots  for  both  modal  and 
zonal  transmitters.  In  Fig.  5-1  the  transmitter  RMS  surface  displacement 
vs.  single  mode  (or  element)  displacement  is  presented.  The  data  is  tabu- 
lated for  comparison  in  Table  5-1.  Note  that  all  nine  MOWACS  modes  have 
steeper  slopes  than  the  zonal  case,  indicating  a stronger  control  over  the 
surface  for  equal  peak  displacement. 

5.1.2  Modal  Vs.  Zonal  Dither  Characteristics 

A fair  comparison  between  the  modal  and  zonal  COAT  systems  requires 
the  same  dither  frequencies  and  system  parameters.  The  difference  between 
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Figure  5-1.  RMS  Phase  Error  of  Transmitter  Surface  for  Both  Modal 
and  Zonal  COAT,  (a)  Modes  1-9.  (b)  7-,  9-,  and 

100-channel  zonal  elements.  Each  mode  (or  element) 

Is  individually  stepped  from  convergence  (phase 
amplitude  ■ <*>peak) » and  surface  statistics  are  cal- 
culated. (Note  different  scale  factor  on  vertical 


TABLE  5-1 

TRANSMITTER  RMS  PHASE  ERROR 


CHANNEL 

Modal 

1 

2 

3 

4 

5 

6 

7 

8 
9 

Zonal 

7 

19 

100 


^RMS  ^ ^peak 


0.678 

0.811 

0.777 

1.14 

0.978 

0.875 

0.892 

0.722 

0.758 

0.594 

0.360 

0.157 
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these  two  systems  is  their  response  to  target  spatial  effects.  As  an  in- 
dividual channel  (or  mode)  is  dithered,  the  resulting  far-field  pattern 
spatially  modulates  the  distributed  radiance  of  the  target  which  propor- 
tionally modulates  the  average  received  signal  in  the  pupil  plane.  The 
results  for  a point  target  are  presented  in  Fig.  5-2  for  comparison  with 
extended  targets. 


Open  loop  results  which  illustrate  the  target  spatial  effects  for 
both  zonal  and  modal  COAT  are  presented  in  Figs.  5-3  and  5-4.  The  normal- 
ized power,  which  is  computed  by  integrating  the  product  of  the  far-field 

pattern  of  the  transmitter  and  the  target  radiance,  is  plotted  versus  the 
mode  amplitude,  or  transmitter  element  phase.  In  all  cases  for  zonal  COAT, 
a decrease  in  the  slope  of  the  curves  (and  hence  dither  signal  level)  occurs 
as  the  target  increases  in  size  and  as  it  becomes  flatter.  Also,  note  that 
the  maximum  power  occurs  at  zero  phase  excursion.  For  the  modal  transmitter, 
there  is  a similar  reduction  in  dither  signal  level;  however,  the  relative 
weighting  of  the  modes  changes  and  for  several  of  the  modes,  the  maximum 
power  does  not  occur  at  zero  mode  amplitude.  The  implication  here  is  that 

equalizing  the  gain,  and  hence  dither  signal  levels,  for  each  of  the  modes 

on  a point  t;  get  does  not  result  in  equal  weighting  for  a distributed 
target.  This  is  a fundamental  difference  between  the  modal  and  zonal  multi- 
dither systems.  In  addition,  depending  on  the  target  brightness  distribu- 
tion, a diffraction-limited  spot  size  on  the  target  may  not  be  obtainable 
for  several  of  the  correction  modes.  The  reason  for  this  is  that  the  trans- 
mitter will  deform  to  maximize  the  received  signal  which  may  not  necessarily 
produce  a diffraction  limited  beam  on  target. 

The  open  loop  results  shown  graphically  in  Figs.  5-3  and  5-4  are 
quantitatively  summarized  in  Tables  5-2  and  5-3.  The  numerical  values  given 
are  the  slopes  of  the  curves  at  a convergence  level  of  957.  (at  convergence 
the  slope  is  zero).  The  amount  of  modulation  of  the  received  signal  for 
each  mode  or  zonal  element  is  directly  proportional  to  the  slope  of  the 
curves  which  depends  on  the  level  of  convergence. 
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a -40  6a 

0 (degrees) 

(a)  Zonal  Transmitter 


100  - ELEMENT  TRANSMITTER 
19  * ELEMENT  TRANSMITTER 

7 - ELEMENT  TRANSMITTER 


m0  i.a- 
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B 

xH  . 

\ 
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4 

1 

0 (Degree*) 

(b)  MOWACS  Transmitter 

Figure  5-2.  Normalized  On  Axis  Power  (Strehl  Ratio)  Vs.  Single  Trans- 
mitter Element  or  Mode  Phase  (from  Convergence).  Over- 
all transmitter  size  is  the  same  for  all  cases. 


<}>  (deg; 

Target  A (dn  ■ 40  cm,  r 


20  cm) 


<J>  (deg) 

Target  B (d_  = 10  cm,  r 


Target  C (dg  = 40  cm,  r£ 


Target  D (dn  = 40  cm,  r =* 


Figure  5-3.  Integrated  Radiance  of  Illuminated  Sphere  Segments  Vs.  Modal 
Amplitude.  Each  of  MOWACS  mode  amplitudes  is  individually 
stepped  in  small  increments  from  0°  to  90°  and  resulting  inte- 
grated power  (far-field  pattern  * target  distribution)  is 
plotted  normalized  to  its  peak  value.  Scenario  parameters  are: 
transmitter  diameter  (dm)  = .6  m,  range  (R)  - 5 km,  X * 10.6  pm 
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TABLE  5-2 

MODAL  TARGET  EFFECTS 


Target 


Point  Target 


Sphere 

Segment 

<1q  = 10  cm 

r - 1 m 
c 


Sphere 

Segment 

dg  * 40  cm 

r = 20  cm 
c 


Sphere 

Segment 

dg  - 40  cm 

r **  1 m 
c 


Sphere 

Segment 

dg  ■ 40  cm 

r * 4 m 
c 


,-5 


NORMALIZED  SLOPE  (10  /deg) 
Mode  No. 


P/P  - 0.95 
max 


*.Mk  (de8) 


1 

173.1 

0. 

2 

371.3 

0. 

3 

332.3 

0. 

4 

517.1 

0. 

5 

460.3 

0. 

6 

400.9 

0. 

7 

416.9 

0. 

8 

332.4 

0. 

9 

345.7 

0. 

1 

167.6 

1.8 

2 

320.8 

9.0 

3 

343.7 

1.8 

4 

469.4 

5.4 

5 

450.0 

0. 

6 

383.8 

0. 

7 

400.5 

1.8 

8 

335.1 

1.8 

9 

329.1 

0. 

1 

149.3 

1.8 

2 

298.2/241.2 

37.8 

3 

295.7 

0. 

4 

389.4 

18.0 

5 

452.5 

0. 

6 

371.1 

0. 

7 

348.4 

0. 

8 

307.2 

0. 

9 

295.3 

0. 

1 

136.6 

1.8 

2 

257.4 

28.8 

3 

296.2 

0. 

4 

354.1 

16.2 

5 

398.2 

0. 

6 

348.8 

0. 

7 

325.8 

0. 

8 

286.3 

0. 

9 

278.1 

* 

0. 

1 

82.2* 

3.6 

2 

95.1 

28.8 

3 

125.7 

0. 

4 

174.2 

0. 

5 

187.4 

0. 

6 

216.7 

0. 

7 

202.4 

0. 

8 

189.8 

0. 

9 

217.9 

0. 

+ 


Indicates  slope  on  both  sides  of  peak; 


Mode  amplitude  extended  greater  than  90" 
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TABLE  5-3 

ZONAL  TARGET  EFFECTS 


Target 


Point  Target 


Number  of 
Channels 


Element 


NORMALIZED  SLOPE  (10  /de; 

p/p  = 0.95 
max  

256.6 

256.6 

148.4 


Sphere  Segment 
d„  = 10  cm 


Sphere  Segment 

dn  = 40  cm 

r = 20  cm 
c 

Sphere  Segment 
dg  = 40  cm 


262.4 

258.9 


268.3 

239.4 


251.5 

213.6 
166.1 
158.2 
118.0 


Sphere  Segment 
d„  = 40  cm 


173.7 

120.7 
126.5 


C - Center  transmitter  element 
M - Middle  ring  element 
0 - Outside  ring  element 


* - Minimum  P/P^ 
t - Minimum  P/P 


0.96  @ <t> 

0.96  @ 4> 
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Dither  modulation  of  the  received  signal  Is  largest  when  the  Illum- 
inating beam  is  unconverged  and  goes  to  zero  at  convergence.  The  amplitude, 
modulation  of  the  received  signal  is  calculated  by  multiplying  the  slopes 
(in  Tables  5-2  and  5-3)  by  the  dither  amplitude  in  degrees.  For  example, 
for  a point  target,  a ±20®  dither  amplitude  for  Mode  1 results  in  the  fol- 
lowing received  signal  modulation:  (slope)  * (dither)/(P/P  ) ■ (173.1  * 

lO"5) (20)/(0. 95)  ■ 0.036,  i.e.,  3.62  modulation  at  95%  convergence. 

5.2  ILLUSTRATIVE  CLOSED-LOOP  SIMULATIONS 

It  must  be  emphasized  that  these  closed-loop  simulations  represent 
illustrative  examples  of  COAT  performance  and  are  not  meant  to  be  a com- 
plete set  from  which  general  conclusions  can  be  drawn.  However,  they  do 
give  an  indication  of  the  performance  with  single-feature  targets  within 
the  limits  of  variation  for  the  statistical  processes  simulated.  They  have 
been  selected  to  reflect  the  trends  in  the  analytic  results  noted  in  Sec.  2. 

The  performance  with  real  targets  depends  on  the  dither  signal  reduc- 
tion by  the  target's  brightness  gradient,  which  is  determined  by  its  geo- 
metrical shape,  size,  and  materials  reflectance,  and  on  the  distribution  of 
the  target  modulation  (noise)  in  the  frequency  domain.  For  zonal  multi- 
dither systems  the  performance  is  reduced  as  the  number  of  transmitter 
elements  (and  hence  control  channels)  is  increased.  Phase  conjugate  per- 
formance may  be  improved  for  more  channels  for  some  target  sizes. ^ 

Modal  multidither  (using  the  first  nine  MOWACS  modes  with  the  same 
interdither  frequency  spacing  and  lowest  dither  frequency)  performance  is 
comparable  to  a 7-element  zonal  system  for  a mode  by  mode  (or  element  by 
element)  Initial  phase  error.  However,  it  achieves  better  performance 
with  the  more  physically  realistic  initial  error  of  an  atmosphere  phase 
screen.  This  is  because  the  modal  shapes  more  closely  approximate  lower 
order  aberrations  than  7 zonal  elements  (in  an  hexagonal  array).  In  fact, 
using  only  the  first  4 MOWACS  modes  provides  correcting  power  superior 
to  that  of  a 7-element  zonal  system  with  respect  to  tilt  and  focus  error. 

\j.  Radley  et  al.,  op.  cit. 
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A 19-element  zonal  transmitter  possesses  approximately  the  same  cor- 
recting power  as  the  9-mode  MOWACS  in  terms  of  atmospheric  phase  perturba- 
tions, but  performs  significantly  worse  with  target  noise  in  the  multidither 
implementation.  It  should  be  noted,  however,  that  the  modal  system  advan- 
tage is  totally  lost  if  higher  order  aberration  correction  is  required.  For 
a dither  system  with  many  channels  (—  50-100)  modal  and  zonal  performance 
should  be  comparable  if  the  dither  frequencies  are  the  same. 

5.2.1  System  Parameters 

There  are  a number  of  design  details  which  must  be  considered  in  order 
to  compare  modal  and  zonal  systems.  For  all  these  closed-loop  runs,  AGC, 
aperture  averaging,  and  disabling  the  error-correction  voltages  to  the 
mirror  pushers  for  the  first  0.125  msec  were  incorporated.  Aperture  averag- 
ing and  AGC  were  found  to  improve  zonal  COAT  performance  by  decreasing  target 
modulation  effects  and  thus  were  incorporated  in  the  modal  simulator  with 
distributed  targets.  It  was  found  that  the  performance  of  the  modal  trans- 
mitter was  improved  by  limiting  the  gain  of  the  AGC  initially  with  a voltage 
clipper.  In  addition,  by  disabling  the  error-correction  voltages  for  the 
first  0.125  msec,  large  initial  transients  in  the  correction  voltages  are 
not  applied  to  the  transmitter.  Otherwise,  the  initial  transients  can 
drive  the  modes  to  a state  which  results  in  a local  maximum  at  convergence. 

Closed-loop  comparisons  between  zonal  and  modal  COAT  require  the  same 
dither  levels.  Therefore,  along  with  the  equivalent  of  the  7-channel  zonal 
dither  level  (±28.8°  peak  amplitude  for  Mode  1 and  the  others  scaled  to  give 
equal  modulation) , an  additional  set  of  runs  was  made  using  the  equivalent 
of  the  19-channel  zonal  dither  level  (±17.1°  for  Mode  1). 

Figures  5-5  and  5-6  illustrate  the  effects  of  target  modulation.  The 
smaller  dither  results  in  a higher  convergence  level  with  a point  target, 
but  performs  worse  with  a distributed  target  when  the  modulation  spectrum 
overlaps  the  dither  frequencies  as  anticipated  from  gain  considerations 
discussed  in  Sec.  2. 
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5.2.2  Modal.  Zonal.  Phase  Conjugate  COAT  Base-Line  Comparison 

As  previously  mentioned,  equivalent  system  parameters  were  applied 
to  the  three  COAT  systems  to  present  a fair  comparison  of  their  performance 
against  varied  target  spatial  effects. 

Typical  results  and  analysis  indicate  that  phase  conjugate  COAT  per- 
forms best  at  fast  rotation  rates,  better  than  multidither  COAT  in  this 
respect.  However,  for  slow  rotation  rates,  multidither  COAT  surpasses 
phase  conjugate.  The  mechanism  for  these  results  is  target  modulation  inter- 
fering with  the  signal  processing.  When  target  base-band  noise  overlaps 
the  dither  bands  for  multidither  COAT,  or  falls  within  the  filter  bandwidth 
of  the  signal  processor  for  phase  conjugate,  erroneous  wavefront  correction 
occurs,  resulting  in  a degradation  in  performance.  Figures  5-7  to  5-8 
illustrate  this  effect.  With  no  target  modulation,  as  in  the  case  of  a 
point  target,  maximum  performance  is  achieved  for  all  COAT  systems. 

Another  very  important  consideration  is  how  the  closed-loop  perform- 
ance scales  with  the  number  of  control  channels.  Other  studies  have  shown 
that  as  many  as  60  to  100  transmitter  elements  (for  a 1-m-diameter  total 
aperture)  may  be  required  to  provide  adequate  compensation  for  atmospheric 
phase  perturbations  along  the  propagation  path.  A modal  transmitter  re- 
quires fewer  control  channels  for  lower  order  aberration  because  each 
dithered  mode  controls  the  entire  transmitter  aperture  which  is  better 
matched  to  the  aberration  than  a zonal  system.  Consequently,  a modal 
transmitter  possesses  higher  correcting  power  for  a small  number  of  channels 
than  an  equal  number  of  zonal  elements.  Therefore,  it  is  not  equivalent  to 
compare  the  9-mode  MOWACS  to  only  the  7-element  zonal  systems.  However, 
when  a large  number  of  channels  is  required  due  to  high  order  aberrations, 
equal  numbers  of  channels  are  essentially  equivalent. 

The  closed-loop  performance  of  a 19-element  zonal  multidither  and 
phase  conjugate  system  is  shown  in  Fig.  5-9.  Note  that  the  multidither 
performance  is  reduced  from  the  7-element  case.  Figure  5-10  shows  the  com- 
parison to  37  elements.  The  reason  for  the  reduction  in  a zonal  multidither 
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Figure  5-8.  Modal,  Zonal,  Phase  Conjugate  COAT  Closed  Loop  Simulations  with  Target  B. 

Curves  are  Strehl  ratio  vs.  time  (0  to  5 msec).  Dots  show  beam  maximum. 


system  is  that  the  dither  modulation  of  the  received  signal  is  reduced  as 
the  number  of  elements  is  increased,  assuming  the  same  overall  transmitter  . 
size  (cf.  Secs.  2 and  3).  Increased  gain  to  recover  adequate  dither 
signals  also  amplifies  the  noise  causing  larger  random  phasing  errors 
In  each  channel.  Phase  conjugate  performance  does  not  suffer  for  larger 
numbers  of  elements  and  in  some  instances  may  actually  improve  as  the 
number  of  channels  increases.* 


The  method  used  to  accomplish  the  19-  and  37-elemunt  zonal  runs  was 
to  generate  the  target  noise  separately  in  an  open-loop  simulation  and  then 
use  it  as  multiplicative  noise  in  the  closed-loop  simulations.  This  is 
valid  provided  that  the  target  is  unresolved  by  the  transmitter  aperture. 
Comparison  to  the  target-in-the-loop  closed-loop  results  and  further  de- 
scription are  given  elsewhere.*  Zonal  and  phase  conjugate  results  are 
shown  in  Figs.  5-11  to  5-13  based  on  analysis  which  indicates  that  trends 
for  a modal  transmitter  will  be  the  same  as  the  zonal. 


The  threshold  process  noted  in  Sec.  2 is  illustrated  clearly  in 
Fig.  5-12  for  both  dither  and  conjugate  systems. 

A 

5.2.3  Closed-Loop  Comparison  with  an  Atmosphere 

Further  performance  comparisons  between  modal,  zonal,  and  phase  con- 
jugate systems  were  made  incorporating  two  static-phase  screen  atmospheres 
of  different  turbulence  levels.  These  systems  were  exercised  against  two 

. 

different  targets:  point  target,  and  a rotating  unresolved  glint  (sphere 

segment,  dQ  = 10  cm,  rc  ■ 1 m).  For  the  distributed  target,  the  most 
stressful  rotation  rate  occurs  when  the  target  modulation  overlaps  the 
dither  band  for  the  two  multidither  systems,  and  for  the  phase  conjugate 
when  TM  falls  within  the  signal  processing  filter  bandwidth.  The  appropri- 
ate targ' t modulation  extents  were  36  and  0.9  kHz,  respectively.  This 
corresponds  to  rotation  rates  of  2 and  0.05  rad/sec.  In  the  modal  COAT 
system,  ±28.8®  (for  Mode  1)  dither  performed  better  than  ±17.7®  dither 
against  target  modulation  interference;  hence  it  was  used  for  these  comparisons. 


-11.  Comparison  of  7-Channel  COAT  Performance  with  Target  C (dg  ■ 40  cm  S 4XR 
Curves  show  Strehl  ratio  vs.  time  (0  to  5 msec).  Dots  show  beam  maximum 


(b)  Phase  Conjugate 


Figure  5-13.  Comparison  of  7-Channel  COAT  Performance  with  Target  A 
(dQ  * 40  era  = 4AR/dT,  rc  = 20  cm).  Curves  show  Strehl 
ratio  vs.  time  (0  to  5 msec).  Dots  show  beam  maximum. 


Without  exception,  the  9-mode  modal  COAT  performed  slightly  better 
than  7-element  zonal  COAT.  This  result  is  expected  since  it  can  closer 
represent  the  ideal  correction  for  a well  behaved  atmospheric  phase  screen 
with  its  deformable  surface. 

Figure  5-14  illustrates  performance  comparisons  for  the  stronger 
atmospheric  turbulence.  For  the  weaker  turbulence  all  three  systems  per- 
formed equivalently  (see  Fig.  5-15),  because  this  phase  perturbation  pos- 
sessed very  slight  tilt  error. 

Closed-loop  performance  runs,  for  a point  target  with  the  phase 
screen  atmosphere,  showed  that  wavefront  corrections  were  mostly  applied 
by  the  first  four  modes  of  the  modal  transmitter,  corresponding  to  focus 
and  tilt.  Illustrative  modal  simulations  indicate  that  much  of  the  atmos- 
pheric perturbation  can  be  corrected  by  these  basic  low  spatial  order  modes 
(see  Fig.  5-16).  It  is  understood,  however,  that  high  spatial  frequency 
atmospheric  perturbations  can  only  be  corrected  using  higher  order  modes  to 
achieve  a high  level  of  convergence. 


The  possibility  of  design  trade-offs  in  scenarios  where  target  noise 
is  a limiting  factor  is  more  substantial  for  modal  than  for  zonal. 

If  the  target  noise  limits  M-D  operation  to  a small  number  of  channels 
C<10) , the  modal  system  would  be  capable  of  providing  greater  turbulence 
correction  and  hence  improved  performance  over  a zonal  system.  Also,  since 
the  modal  system  involves  whole  aperture  control,  an  adaptive  system  which 
tailors  the  number  of  channels  and  gain  parameters  to  optimize  power  on 
target  could  be  more  easily  implemented  than  a zonal  system.  However,  these 
options  are  of  interest  only  if  optimization  at  less  than  diffraction 
limited  Intensity  levels  is  adequate  to  satify  the  mission  requirements. 
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time  (0  to  5 msec). 


APPENDIX  A 

MODAL  TRANSMITTER  SIMULATION 

In  this  section  Is  described  a general  mathematical  formalism  for 
representing  the  vibrant  structural  modes  of  a deformable  mirror  and  how 
this  is  implemented  on  a digital  computer  to  simulate  a modal  adaptive 
transmitter.  The  particular  choice  for  this  study  is  the  MOWACS  (MOdal 
Wavefront  Adaptive  £ontrol  System)  mirror  developed  by  Perkin-Elmer 

In  the  first  part  of  this  appendix  we  present  the  method  of  describ- 
ing the  first  nine  MOWACS  modes  in  terms  of  a series  of  Zemike  polynomials. 
In  the  second  part  some  of  the  practical  considerations  are  given  for  rep- 
resenting the  transmitter  for  both  the  open-  and  closed-loop  computer  sim- 
ulations presented  in  Sec.  5. 

A.l  METHOD  OF  MODELING  MOWACS  MIRROR 

A deformable  mirror  which  is  driven  in  one  or  more  resonant  struc- 
tural modes  by  an  array  of  actuators  is  called  a modal  mirror.  For  prac- 
tical reasons,  this  type  of  mirror  is  usually  clamped  at  two  opposing  edges 
to  accommodate  coolant;  the  other  two  edges  are  free.  The  various  modes  of 
vibration  (usually  sinusoidal)  that  can  be  sustained  are  excited  by  pushers 
located  on  the  underside  of  the  mirror.  Virtually  any  set  of  orthogonal 
polynomials  can  be  used  to  describe  the  mirror-surface  displacement  of  each 
mode  and  can  be  summed  to  give  the  total  displacement  as  a function  of  time 
and  position  across  the  surface. 

The  mathematical  representation  used  by  Perkin-Elmer  for  the  spatial 
dependence  of  the  mirror  surface  displacement  is  the  first  22  Zemike  poly- 
nomials (P j ) which  are  listed  in  Table  A-l.  The  coefficients,  a^ , that 
accompany  these  polynomials,  in  the  description  of  the  first  nine  modes  of 
vibration  are  listed  in  Table  A-2.  Thus,  the  mirror  surface  displacement 
h^(x,  y)  of  the  ith  mode  of  vibration  is  approximated  as 

^C.  Neufeld,  Modal  Wavefront  Control  System  (MOWACS)  Phase  II  Interim 
Report,  Perkin-Elmer  Corporation  Report  12700,  September  1975. 
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MODIFIED  ZERNIKE  POLYNOMIALS* 


P 2 = 2rcosQ  = 2x 
» 2rsin0  = 2y 
P4  = y3(2r2-l) 

I 2 

Pj  « f6  r sin20  = y6  2xy 
p6  * v6  r2cos29  - v6  (x2-  y2) 

3 2 

P?  = y8  (3r  - 2r)  sin0  = '/8.  y(3r  - 2) 

P = V"8  (3r3-  2r)  cos0  = \f8  «x(3r2-  2) 

P9  = y8  r3  sin3Q  = /8  y(3x2-  y2) 

P1Q  = V"8  r3  cos39=  y8  x(x2-  3y2) 

Pu  = V5  (6r4-  6r2+  1) 

pi2 = 1,ri°  (4r4* 3r2)  cos20  = V"10  (x2‘  y2H4r2- 3) 

P13  = yio  (4r4-  3r2)  sin29  = yiO  2xy  (4r2-  3) 

P^4  = yiO  r4  cos4Q  = yiO  (x4-  6x2y“+  y4) 

P^3  = yiO  r4  sin4Q  = y 10  4xy  (x2-  y2) 

P16  = yi2  (10r5-  12r3+  3r)  cosQ  = yi2  x(10r4-  12r2+  3) 
P1?  = yi2  (10r5-  12r3  + 3r)  sin9  = yi2  y(10r4-  12r2+  3) 
P18  = yi2  (5r5-  4r3)  cos39  = /12  x(x2-  3y2)(5r2-  4) 

P19  = yi2  (5r5-  4r3)  sin30  = y 12  y.(3x2-  y2)(5r2-  4) 

P2Q  * yi2  r5  cos50)  = yi2  x(x4-  10x2y2+  5y4) 

P21  = V 12  rJ  sin50  = yi2  y(5x4-  10x2y2+  y4) 

P22  = (20r6_  30r4+  12r2"  1) 

* 

Reproduced  from  Neufeld,  op.  cit. 
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TABLE  A-2 

ZERNIKE  DECOMPOSITION  OF  9 MODE  EXPERIMENT  MODES’ 


L 

j 

j 


\ 


CM 


O 

H 

cm 


ON 

CM 


QO 

CM 


vO 

tM 


m 

CM 


CM 


m 

CM 


CM 


CM 


CM 


O 

•o  • 

o o 


* 

o 

in 

o 

© 

O 

o 

© 

r> 

o 

93 

o 

© 

o 

© 

mm 

o 

o 

IT 

© 

o 

o 

© ® 

© 

o 

o 

o 

© 

© r» 

© 

mm 

© 

O 

o 

© 

o 

o 

3 

© 

© 

o 

© 

o 

o 

© 

o 

•• 

o 

o o o 

© © © o © © 

1 

« 

f 

• 

o 

o 

c 

o 

© 

mm 

MOO 

o 

c 

o 

Aj  © 

C 

o 

O 

o 

o 

o © 

mm 

o 

•9 

o 

o o 

o 

o 

o 

9 

© 

o 

© 

© 

o 

o c 

X 

o 

-c 

o 

o 

© 

© 

o 

o 

mm 

o 

AJ 

© 

© 

© 

o © o 

o 

© 

© 

o © 

o 

• 

• 

1 

• 

1 

• 

• 

c rvj  c 

© aj 

O 

o 

o 

© e 

o 

— • 

o 

© 

© 

o o o 

O O' 

© 

© 

© 

o 

p* 

o 

© 

© 

o 

© 

o a- 

© 

o rN 

© 

© 

© 

o 

3 

© 

© 

3 

© 

© AJ 

o 

© 

© 

o 

o 

O 

o 

3 O 

1 

1 

1 

1 

1 

c 

c 

c 

© X 

o 

© 

© 

© o 

c 

<n 

© © 

o © 

© 

© 

c 

© 

© 

© 

© 

c 

© 

© 

o 

o 

© 

** 

© 

o 

© 

o 

© 

3 

o 

© aj 

O 7 

o 

© 

o 

© 

o 

© 

O Aj 

© 

© 

© 

o o o 

© 

O 

o 

© o 

o 

1 

1 

1 

1 

1 

o 

© 

© 

o 

© 

© 

<n 

O K 

© 

© 

Aj 

o 

c 

o 

■9 

o 

p— 

c 

C 

r 

c 

© 

© *■“ 

o 

c 

© 

o 

© 

© 

c 

m 

© 

O Aj 

o 

o 

© 

Aj 

© 

c 

c 

© r> 

© 

3 

© 

© 

o 

© 

c 

© 

3 © 

0 

0 

0 

• 

1 

1 

1 

1 

in  o a 

© A| 

o 

o 

© 

<n  o 

9 

© 

3 © 

© 

-r 

3 

* 

c 

c 

C 

© 

© 

© 

© 

© 

9 

© r- 

© 

3 

o 

3 

O 

c 

9 

c 

3 

3 

© 

n 

3 

*• 

O Aj 

o 

O 

© 

3 

3 

O 

o o 

© 

0 

0 

0 

© 

O 

O 

1 

• 

f 

1 

t 

1 

o o © 

© 

© 

AJ  © 

mm 

o © 

o 

o 

© 

3 

* 

O 

O 

o c 

© 

© 

© 

© 

o 

© 

© 

o 

o 

o 

© 

3 

3 

n © 

© 

© 

o c 

© 

O 

© 

o 

O 

c 

© 

c 

o 

© 

© 

o 

c 

c 

© 

3 

o 

© 

o o 

© 

• 

1 

• 

1 

• 

< 

© 

© 

© 

© 

o 

<r 

o 

< 

© 

7 

© 

e 

O 

o 

3 

o 

— 

c 

9 

C 

© 

o 

o 

x 

© 

< 

© 

x 

© 

© 

© 

c 

•4 

© 

© 

•9 

© 

c 

c 

3 

r> 

© 

o m 

© 

© 

© 

o 

r-> 

o 

c 

c 

C 

© 

o 

© 

o 

3 

© 

I 

« 

» 

1 

1 

AJ  «T 

p«* 

© c* 

3 

© 

© 

© 

P-. 

C 

© 

-n 

C 

© 

© X 

e 

n 

© 

© 

•m 

C 

o 

3 

mm 

o 

3 

c 

o 

>9 

3 

© 

3 

9 

c 

3 

o 

o 

© 

© 

9 

o 

c 

c 

- 9 

© 

© 

© 

3 

e c 
1 

o 

© 

3 

o 

© 

O 

e 

© 

X 

o «C 

j , 

© 

c 

c 

e 

C7* 

© in 

c 

O 

© 

c 

e 

© 

© 

3 

o 

o 

© 

e 

c 

c 

p* 

c 

mm 

c 

o 

© 

c 

© 

© 

IT 

© 

mm 

c 

3 

O 

c 

c 

© 

•T 

c 

r> 

o 

c 

3 

c 

c 

© 

C 

© 

- 

c 

© 

© 

l 

1 

1 

1 

© 

7 

© 

© 

© 

c 

© 

ru 

c 

in 

© 

© 

o 

•3 

mm 

3 

o 

c 

x 

C 

c 

c 

c 

X 

O 

9 

c 

n 

© 

•<5 

© 

© 

x 

© 

r' 

c 

c 

o 

c 

o 

n 

© 

© 

m 

3 

© 

© 

c 

© 

o 

• 

• 

• 

• 

• 

* 

• 

• 

• 

c 

c 

© 

C 

© 

c 

© 

c 

c 

1 

1 

i 

r* 

-0 

iT 

< 

K 

«* 

7 

CM 

CM 

CM 


O 

CM 

CM 


O 

CM 


00 

CM 


m 

pH 

CM 


cn 

r-H 

CM 


CM 

H 

CM 


in 

© 

© 

o o 

o 

A- 

o 

4 

O 

7 

© 

o 

o 

o 

P4 

o 

n 

O 

m 

o 

o 

o 

O AJ 

o 

o 

o 

PM 

o 

o 

o 

© <\J 

© 

© 

o 

© 

© 

o 

© 

O Aj 

o 

e 

© 

o 

o 

o 

o 

o 

© 

© 

© © © 

o 

o 

o 

0 

0 

0 

• 

1 

1 

1 

« 

« 

C 

© 

© 

o 

mm 

© 

o 

o n 

o 

<* 

o 

o 

mm 

o 

© 

o in 

© in 

o 

© lT 

o 

o 

© 

mm 

©n© 

© 

mm 

o 

o 

o 

9 

o 

© 

© 

O CSJ 

© 

o 

O p* 

o 

o o 

o 

© 

© 

o 

o 

o 

© 

o 

o 

o 

o 

o 

© 

© 

1 

• 

1 

1 

• 

1 

© © 

© 

3 

© 

© © © 

© 

© 

© 

p-a 

o 

> 

o 

© 

O 

© 

o 

© 

© 

© *M 

© 

o 

o 

© 

o 

3 

© 

o 

o 

© 

© 

© 

O 

mm 

© ^ 

o 

© 

© 

3 

© 

O 

© 

© 

o 

o 

o 

© 

0. 

0 

o 

© 

© 

© 

© 

© 

o 

1 

• 

• 

1 

• 

1 

O 

7 

c 

© 

9 

© 

© 

o O 

© in 

© 

© 

mm 

3 

3 

o 

© 

© 

3 

c 

© 

T 

3 

o 

© 

© 

© 

© 

X) 

O 

© 

© 

© 

© 

© 

© 

mm 

o 

© 

O AJ 

c 

© 

© 

© 

© 

o 

© 

o 

© 

© 

o 

o 

© 

o 

<3 

© 

o 

© 

• 

• 

1 

1 

< 

1 

• 

o 

o 

© 

9 

© 

3 

mm 

© 

o 

c 

© 

© 

7 

© CVI 

c 

© 

© 

© 

© 

© 

9 

O M 

© 

o 

© 

o 

c 

s 

f\j 

3 

O 

© 

© 

© 

© 

c 

c 

M 

© 

p*» 

o 

o 

© 

3 

© 

c 

© 

3 

© 

o 

© 

o 

O 

© 

o 

3 

© 

© 

o 

© 

© 

1 

1 

1 

m 

© 

X 

© 

o 

«« 

p* 

3 

3 

3 

c 

© 

3 

3 

AJ 

3 

\ 

C 

© 

■n 

© 

© 

3 

X 

s 

\i 

© 

c 

3 

© 

3 

3 

C 

© 

3 

3 

O 

3 

© 

© 

© 

c 

© 

O 

O 

o 

O 

• 

• 

• 

• 

• 

• 

• 

• 

• 

3 

o 

© 

© 

O 

© 

3 

3 

© 

« 

« 

1 

• 

1 

O 

© 

© 

_ 

o 

iT 

© 

© 

© 

3 

3 

c 

© 

3 

© 

3 

3 

3 

O 

© 

X 

o 

r-» 

3 

© 

© 

3 

o 

3 

•m 

3 

O 

3 

O 

O 

3 

© 

C 

<M 

3 

9 

O 

© 

3 

O 

© 

© 

C 

3 

O 

© 

3 

© 

• 

• 

• 

• 

• 

• 

• 

• 

• 

© 

© 

© 

© 

© 

o 

o 

o 

O 

1 

1 

1 

1 

1 

1 

o 

3 

© 

-i 

C 

c 

e 

© 

© 

© 

© 

o 

3 

3 

3 

O 

o 

© 

C 

© 

C* 

c 

© 

X 

3 

3 

3 

3 

© 

© 

© 

9 

C 

3 

o 

© 

c* 

3 

3 

7 

© 

© 

c 

O 

© 

© 

3 

3 

© 

3 

O 

© 

e 

© 

c 

© 

© 

3 

O 

© 

• 

• 

1 

1 

• 

<r 

c 

3 

© 

© 

3 

© 

— 

aj 

c 

© 

© 

© 

© 

© 

r- 

© 

© 

C 

3 

© 

r- 

3 

rm 

o 

v» 

© 

rr 

3 

3 

© 

© 

mm 

© 

T 

© 

o 

O 

O 

X 

O 

© 

© 

© 

© 

C 

© 

o 

© 

© 

0. 

n. 

c 

© 

© 

o 

© 

o 

1 

1 

1 

1 

© 

© 

C 

© 

© 

C 

* 

© 

e 

c 

C 

C 

© 

© 

C 

in 

o 

© 

C 

— 

© 

Aj 

© 

© 

C 

© 

c 

c 

~ 

3 

* 

© 

© 

© 

© 

© 

c 

© 

O 

X 

o 

o 

© 

© 

© 

o 

© 

C 

A. 

© 

© 

© 

c 

© 

e 

C 

© 

© 

© 

o 

• 

• 

1 

1 

1 

1 

© 

© 

X 

c 

© 

c 

3 

AJ 

o 

r* 

© 

mm 

© 

o 

3 

© 

o 

X 

c 

c 

c 

C 

3 

r* 

o 

n 

© 

ru 

© 

c 

© 

o 

A. 

o 

mm 

c 

_ 

c; 

O 

© 

p-» 

o 

© 

© 

c 

© 

© 

c 

c 

C 

© 

• 

• 

• 

• 

• 

• 

• 

• 

• 

© 

© 

c 

1 

c 

I 

c 

c 

c 

© 

© 

mm 

n. 

r-i 

w 

X 

© 

72 


Reproduced  from  Neufeld,  op.  cit.  The  Z,  coefficients  for  the  i,th  mode  are  the  in  Eq.  A-l 


h. (x,  y)  - £ 
* 


a P 

ij  j 


(A-l) 


The  total  mirror  surface  displacement  Is  the  sum  of  the  h^,.  The  approxi- 
mation of  the  surface  displacement  by  Zernike  polynomials  is  quite  accurate 
provided  (x,  y)  lies  within  a unit  circle,  the  range  over  which  the  Zernike 
polynomials  are  defined.  Accordingly,  the  mirror  is  mathematically  masked 
to  exclude  these  non-defined  regions.  The  maximum  amplitude  of  the  dis- 
placement of  any  mode  is  normalized  to  unity. 

This  Zernike  polynomial  formalism  is  conveniently  represented  as  a 
transformation  matrix  for  incorporation  into  the  existing  COAT  computer 
simulators.  This  matrix,  called  the  error  signal  processor  (cf.  Fig.  A-l), 
calculates  the  total  mirror  displacement  h(x,  y)  for  any  combination  of  indi- 
vidual mode  amplitudes.  The  processor  has  nine  columns  corresponding  to  each 
of  the  MOWACS  modes  and  n rows  corresponding  to  the  number  of  points  describ- 
ing the  mirror  figure. 

Sample  results  generated  by  the  error-signal  processor  are  shown 
graphically  in  Figs,  A-2  and  A-3.  In  the  first,  each  of  the  mode  shapes 
is  displayed,  ana  in  the  second,  the  time  evolution  of  modes  4 and  5 is 
shown. 


The  Zernike  polynomials  used  in  describing  the  mirror  surface  dis- 
placement were  chosen  not  for  mathematical  simplicity,  but  because  various 
ones  are  associated  with  aberrations  of  the  resulting  optical  wavefront  of 
the  reflected  beam  (Fig.  A-4) . The  expression  for  each  aberration  shown  is 
for  a Zernike  polynomial  given  in  Table  A-l  and  differ  somewhat  from  the 
definition  of  the  classic  aberrations.  The  reason  is  that  an  additional  set 
of  aberrations  that  is  equivalent  to  the  set  given  here  exists,  but  they 
have  a rotated  axis  of  symmetry  whose  rotation  depends  on  the  order  of  the 
aberration.^-  This  is  merely  a choice  of  nomenclature  and  does  not  affect 
the  results  presented. 

Hi.  Born  and  E.  Wolf,  Principles  of  Optics,  Pergamon  Press  Ltd.,  1970, 

p.  212. 
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;ure  A- 2.  First  Nine  MOWACS  Modes.  Plots  illustrate  magnified  displacement  only  within  circular 
illuminated  area. 


Evolution  of  Superposition  of  Modes  4 and  5 (Frequencies  6214  Hz,  7261  Hz,  Respec 
y)  In  Intervals  of  18°  Phase  of  Mode  4. 


A.  2 SIMULATOR  DESIGN  FOR  TRANSMITTER 

Unlike  the  zonal  transmitter,  the  modal  transmitter  has  a continuous  . 
deformable  surface  which  requires  more  transmitter  plane  points  to  repre- 
sent the  modal  surface  figure  than  the  zonal  case.  It  would  be  to  our 
advantage  to  minimize  the  number  of  transmitter  point  locations  and  yet 
accurately  model  the  modal  surface  from  the  computer  simulator  run  time 
point  of  view. 

A. 2.1  Sampling  Theorem 

Although  the  MOWACS  vibrational  modes  are  a summation  of  Zernike 
polynomials  which  have  different  spatial  frequencies,  this  may  not  deter- 
mine the  highest  overall  modal  spatial  frequency.  Observing  Fig.  A-2,  it 
appears  that,  in  terms  of  trigometric  functions,  the  largest  number  of 
cycles  across  the  transmitter  is  not  more  than  two.  Nyquist  theorem  states 
that  the  minimum  sampling  frequency  of  any  signal  should  be  at  least  two 
times  the  highest  frequency  component  of  that  signal. 

A. 2. 2 16  x 16  7s.  8x8  Transmitter 

Because  of  the  edge  effects  of  the  transmitter  and  the  uncertainty  of 
accurately  modeling  them,  two  transmitter  configurations  were  considered 
with  larger  grids  (16  x 16  and  8 x 8)  than  required  by  the  Nyquist  cri- 
terion. 

Vibrational  modes  were  compared  in  the  far-field  target  plane  to 
determine  the  degradation  of  the  8 x 8 transmitter  over  the  16  x 16.  In 
Fig.  A-5,  Mode  7 was  applied  to  both  transmitters  and  a comparison  was  done 
in  the  far  field.  For  a tt/2  rad  modal  amplitude,  the  standard  deviation 
(RMS)  of  the  error  observed  was  less  than  17.  with  a maximum  distortion  of 
5.7%  (0.046  + 0.011).  Applying  all  modes  with  random  modal  amplitudes 
ranging  from  -rr/2  to  tt/2  produces  an  RMS  error  of  1.7%  and  a maximum  dis- 
tortion of  13%  (cf.  Fig.  A-6) . Ba^ed  on  these  results  the  8 x 8 grid  was 
considered  acceptable. 
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Figure  A-5.  Error  Calculation  for  Representing  MOWACS  Mode  //7  by  64  Vs.  256  Points.  Transmitter 
phase  error  is  n/2.  In  = maximum  (diffraction  limited)  intensity. 


A.2.3  Target  Plane  Beam  Patterns 

A series  of  target  plane  beam  profiles  (shown  in  Figs.  A-7  through 
A-9)  were  calculated  for  individual  and  combinations  of  modes  for  various 
amplitudes  to  determine  the  distribution  of  power  on  the  target.  Because 
each  mode  controls  the  entire  transmitter  aperture,  dither  amplitudes 
equal  to  those  used  in  zonal  systems  give  significantly  larger  fluctua- 
tions in  target  plane  power  distribution  than  a zonal  transmitter. 


Lgure  A-7.  Target  Plane  Beam  Patterns  for  MOWACS  Mode  #1.  Vertical  displacement  is  In- 
tensity (linear  scale)  normalized  to  diffraction  limited  (<p  = 0) . Scenario 

parameters  are:  range  = 5 km,  X = 10.6  pm,  (J>  =*  transmitter  phase  error. 
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APPENDIX  B 

MODAL  BASE-LINE  SYSTEM  EVALUATION 


B.l  MODAL  CORRECTION  MECHANISM 

Although  initial  perturbations  are  applied  to  the  transmitter  wave- 
front  mode  by  mode,  it  is  difficult  to  determine  the  mechanism  by  which 
the  overall  modal  perturbation  is  corrected.  It  appears  that  a combina- 
tion of  both  mode  by  mode  and  summation  of  modes  approximate  the  conjugate 
of  the  initial  perturbation.  In  other  words,  each  initial  mode  perturba- 
tion is  corrected  by  the  negative  of  that  perturbed  mode.  In  addition,  a 
combination  of  correcting  modes  may  approximate  the  necessary  correction 
to  the  wavefront,  but  not  necessarily  mode  for  mode.  The  overall  effect 
appears  to  be  the  summation  of  all  correcting  modes  working  together  to 
correct  for  the  overall  initial  phase  perturbation.  Figure  B-l  illustrates 
this  effect,  where  the  resultant  phase  error  (initial  mode  amplitude  plus 
correcting  mode  amplitude)  is  plotted.  In  case  (a)  an  initial  phase  per- 
turbation is  applied  only  to  Mode  9.  Note  that  all  modes  are  correcting 
for  the  perturbation  caused  by  Mode  9.  In  case  (b) , both  Modes  4 and  9 
are  perturbed  initially.  In  both  cases,  the  mode(s)  that  was  initially 
perturbed  contributed  the  maximum  correction,  but  additional  correction 
is  provided  by  the  other  modes. 

B.2  SENSITIVITY  TO  INITIAL  TRANSMITTER  PHASE  ERROR 

The  additive  property  of  the  different  modes  causes  an  initial 
randomization  of  the  modes  that  can  produce  a wavefront  distortion  too 
severe  for  the  system  to  correct.  Hence,  proper  treatment  in  determin- 
ing the  initial  wavefront  distortion  must  be  exercised.  This  effect  is 
demonstrated  in  Fig.  B-2,  in  which  the  two  target  plane  beam  profiles 
were  produced  by  initially  unconverged  transmitters  with  the  same  RMS 
(1.7  rad)  but  slightly  different  peak-to-peak  errors  (6.28  rad  vs.  6.50 
rad).  The  closed-loop  performance  against  a point  target  (Fig.  B-3a) 
shows  the  sensitivity  of  the  modal  performance  to  the  assumed  initial 
phase  error. 
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Figure  B-l.  Mod3l  Multidither  COAT  Closed  Loop  Performance  with  a Point 
Target  for  One-  and  Two-Mode  Initial  Errors.  Curves  show 
Strehl  ratio  (I/It))  and  amplitude  error  per  mode  (<pei)  vs. 
time.  Modal  transmitter  consists  of  the  first  9 MOWACS 


i Profiles  for  Two  MOWACS  Transmitters  with  Equivalent  Phase  RMS 
Different  Peak-to-Peak  Wavefront  Distortion  (4>p_p) . Parameters  are 
10.6  pm,  transmitter  diameter  = 0.6  m. 
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These  same  Initial  phase  errors  show  even  more  pronounced  effect 
on  a distributed  target  (cf.  Fig.  B-3b) . The  target  is  Target  B (a  sphere 
segment  of  10-cm  extent,  1-m  radius  of  curvature,  comprised  of  aircraft 
material)  rotating  at  0.25  rad/sec.  The  target  size  corresponds  to  half 
the  converged  beam  diameter.  The  bandwidth  of  the  target  noise  is  4.5  kHz, 
which  is  well  below  lowest  dither  frequency  (8  kHz) . The  results  with  the 
more  favorable  initial  beam  profile  show  a similar  level  of  performance 
to  the  zonal  multidither  system. 

B.3  EFFECT  OF  DITHER  FREQUENCIES  ON  SYSTEM  PERFORMANCE  WITH  A DISTRIBUTED 

TARGET 

A fair  comparison  for  both  multidither  systems  require  the  same  dither 
frequencies  and  dither  frequency  separation  if  the  same  electronic  signal 
processor  is  used.  As  dither  frequencies  are  lowered,  base-band  target 
noise  overlaps  the  dither  bands  for  slower  target  rotations.  The  result 
is  a degradation  in  performance  due  to  target  noise  overlapping  the  lower 
dither  frequencies.  In  addition,  if  the  dither  frequency  spacing  is 
decreased,  it  can  be  shown  that  the  convergence  will  be  delayed. 

In  Fig.  B-4,  a comparison  in  modal  performance  for  two  dither  bands, 

the  zonal  equivalent  (8-20  kHz)  and  the  actual  MOWACS  frequencies  (2.3-13.3 

kHz)  is  presented.  The  target  is  Target  B rotating  at  1 rad/sec  (f  = 

max 

18  kHz).  As  expected,  qualitative  results  indicate  a degradation  in  per- 
formance, both  in  convergence  level  and  convergence  time. 

With  the  point  target,  a delay  in  convergence  time  was  expected  for 
the  actual  MOWACS  dither  frequencies  (2.3-13.2  kHz)  due  to  the  additional 
intermodulation  terms.  In  practice,  this  was  found  not  to  be  significant 
as  illustrated  in  Fig.  8-5,  where  similar  performance  between  equally 
spaced  dithers  and  actual  MOWACS  resonant  frequencies  is  compared. 
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(b)  = 2.3-13.3  kHz  (actual  MOWACS  frequencies) 


Figure  B-4.  Modal  Closed  Loop  Performance  with  Target  B for  Two 
Different  Sets  of  Dither  Frequencies.  In  (a)  the 
dithers  are  equally  spaced  Af  apart.  Parameters: 
number  of  modes  =>  9,  range  = 5 km,  \ = 10.6  ym,  = 
1 rad/sec  (f  = IS  kHz). 
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APPENDIX  C 

ATMOSPHERIC  TURBULENCE  MODEL 


The  effect  of  atmospheric  turbulence  on  a received  laser  beam  wave- 
front  is  included  in  system  calculations  by  generating  a two-dimensional 
array  of  phase  fluctuations  induced  by  the  atmosphere  at  the  transceiver 
plane.  The  inclusion  of  such  an  array  here  accounts  for  most  of  the  phase 
variations  caused  by  propagation  through  the  turbulent  atmosphere.  While 
not  including  amplitude  variations  due  to  propagation  from  turbulent  layers 
at  high  altitudes,  such  effects  are  of  secondary  importance  for  the  laser 
system  and  scenarios  of  interest. 

The  array  of  phase  fluctuations  is  computed  from  the  convolution  of 
a set  of  zero  mean,  normally  distributed  random  numbers  with  the  autocorre- 
lation of  the  spatial  frequency  distribution  for  the  propagation  path  under 
consideration.  This  autocorrelation  is  the  Fourier  transform  of  the  spatial 

frequency  spectrum  (assumed  to  be  modified  von  Karman) . For  the  simulation, 

-4 

the  spectrum  was  approximated  by  a < dependence  (where  < is  the  spatial 
frequency)  with  a sharp  cutoff  at  < = corresponding  to  the  outer 

scale  of  turbulence.  The  variance  for  the  phase  fluctuations  derived 
below  is  given  by 

o2  **  0.0087k2C2  dL„5/3  (C-l) 

n n u 

where  k ■ wave  number  2tt/X 
2 

Cn  * refractive  index  structure  parameter 
Lq  ■ outer  scale  of  turbulence 
d ■ effective  width  of  phase-perturbing  region 


^J.  Gilbert  et  al.,  Sensor  Modulation  Effects,  General  Research  Corpora- 
tion, 905-01-CR,  to  be  published. 
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2 

The  square  root  of  Eq.  C-l  is  plotted  in  Fig.  C-l  vs.  Cn  for  sev- 
eral values  of  d and  Lq.  As  is  illustrated  by  the  curves,  a particular 
value  for  the  phase  variance  can  be  obtained  from  numerous  sets  of  the 
other  three  parameters. 


A great  deal  of  experimental  work  has  been  performed  to  determine 
2 

the  parameter  C . For  reference.  Fig.  C-2  presents  representative  experi- 
1 n 2 

mental  data  on  the  variation  of  C as  obtained  from  thermosonde  measure- 

n 

ments . 


The  expression  for  0 follows.  The  refractive  index  spectrum  hy- 
2 n 

pothesized  by  Tatarski  has  the  following  form: 

4>n(<)  - 0.033C2k11/3  exp  (-<2/k2)  (C-2) 

2 

where  C ■ index  of  refraction  structure  function 
a 

K m spatial  frequency 

< m 5. 92/ in  which  ■ inner  scale  of  the  atmosphere 
mu  u 

In  the  limit  as  < approaches  0,  this  spectrum  becomes  indeterminate. 

3 

Strohbehn  added  an  additional  term  that  causes  the  spectrum  to  saturate 
for  values  of  < < 1/Lq.  This  modified  spectrum  has  the  form 

4*1  (k)  - 0.033C2 [k2  + 4n2L  “2]"11/6  exp  (-k2/k2  ) (C-3) 

n n U \ m / 


r.  0.  Mlnott  et  al..  Results  of  Balloon  Atmospheric  Propagation  Experi- 
ment Flights  of  1970.  Goddard  Space  Flight  Center  Document  X-524-72-125 
(March  1972). 

2 

V.  I.  Tatarski,  The  Effects  of  the  Turbulent  Atmosphere  on  Wave  Propaga- 
tion, Report  TT-68-50464,  U.S.  National  Technical  Information  Service, 
Springfield,  Virginia,  p.  76. 

3J.  W.  Strohbehn,  Proc.  IEEE.  56,  p.  1301  (1968). 
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(rad) 


Figure  C-l. 


Curves  of  RMS  Phase  (an)  of  Fluctuations  vs.  the  Refractive 
Index  Structure  Parameter  (C„)  Showing  the  Effects  of  Model 
Parameters  on  the  Magnitude  of  the  Phase  Variations:  Lg  * 

outer  scale  of  turbulence,  d * effective  width  of  phase 
screen 


. j 
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Figure  C-2.  Empirical  Modeling  of  the  Atmosphere 


As  demonstrated  by  Taylor,^-  the  variance  of  the  index  of  refraction  varia- 
tions is  obtained  by  the  following  equation: 


2 2 
ir  k d 


00 

f *(<) 

n 


K d< 


(C-4) 


where  d = width  of  the  phase  screen.  This  integral  is  tabulated,  from 
which  we  obtain,  for  the  condition  Lq  » £q,  Eq.  C-l: 


2 2 5/3 

0.0087k dC  Ln ' 
n 0 


S.  Taylor  et  al.,  JOSA,  65,  p.  78  (1975). 
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APPENDIX  D 

EFFECT  OF  APERTURE  AVERAGING  ON  SPECKLE  POWER  SPECTRAL  DENSITY 

COAT  systems  are  degraded  by  target  effects  because  the  random  phase 
and  intensity  pattern  that  is  generated  by  the  microscopically  irregular 
target  surface  manifests  itself  as  a noise  signal  at  the  COAT  receiving 
aperture.  Accordingly,  it  is  important  to  know  the  extent  and  distribution 
in  the  frequency  domain  of  this  target  modulation  signal  in  comparison  with 
the  COAT  sensing  and  bandpass  frequencies. 

There  are  two  methods  available  to  assess  the  COAT  degradation  by 
target  noise.  The  direct  method  is  the  more  laborious:  the  spatial  distri- 

bution of  target  brightness  is  represented  by  a system  of  point  scatterers 
whose  speckle  pattern  (intensity)  at  the  receiving  plane  is  calculated 
point  by  point  in  its  propagation  from  target  to  receiver;  the  square  of 
the  absolute  magnitude  of  the  Fourier  transform  of  the  speckle  pattern 
gives  the  power  spectral  density;  and  a large  number  of  samples  provides 
the  ensemble-averaged  power  spectral  density  of  the  speckle  pattern  at  the 
receiver. 

The  second  method  to  assess  target-induced  noise  effects  is  analyt- 
ical and  circumvents  a large  amount  of  computation  that  is  required  by  the 
direct  method.  It  is  simply  that,  by  using  properties  of  the  Fourier 
transforms  that  are  involved,  the  power  spectral  density  of  the  intensity 
in  the  receiving-aperture  plane  is  proportional  to  the  autocorrelation  of 
the  target  brightness  distribution.  The  only  assumptions  necessary  are  that 

the  target  and  receiver  be  in  Fourier  transform  planes  and  that  the  target 

12  3 4 

be  diffuse.  Several  authors  * * ’ have  derived  this  relationship  between 
speckle  power  spectral  density  (PSD)  and  target  brightness,  so  only  the  re- 
sults will  be  summarized  here. 

1R.  I.  Goldf ischer , JOSA,  55,  247  (1965). 

2J.  W.  Goodman,  Proc.  IEEE.  53,  1688,  (1965). 

3R.  B.  Crane,  JOSA,  60,  1658  (1970). 

AP.  H.  Dietz,  JOSA,  65,  279  (1975). 
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The  first  assumption  is  that  the  target  and  receiver  electric  fields 
are  related  by  Fourier  transform.  This  is  valid  if  either  one  of  two  con- 
ditions is  met:  (1)  a lens  is  present  in  the  receiver  plane  whose  focal 

length  is  equal  to  the  target  range,  or  (2)  the  target  is  in  the  far  field. 

In  this  report,  perfect  autofocus  is  assumed  and  condition  IH  is  satisfied. 
One  method  of  deriving  the  desired  relationship  is  using  the  autocorrelation 
(or  Wiener-Khintchine)  theorem  of  Fourier  transform  theory. ^ Displayed 
graphically  in  Fig.  D-l,  it  states  that  the  power  spectral  density  of  an 
arbitrary  function  may  be  calculated  either  by  (1)  the  square  of  the  abso- 
lute magnitude  of  the  Fourier  transform  of  the  function,  or  (2)  the  Fourier 
transform  of  the  autocorrelation  of  the  function.  Repeated  application  of 
the  autocorrelation  theorem  to  the  electric  field  distribution  at  the  target 
results  in  the  assorted  relationships  that  are  displayed  in  Fig.  D-2. 

The  final  step  in  relating  the  speckle  PSD  to  the  target  brightness 

distribution  is  to  derive  an  expression  for  the  intensity  autocorrelation 

in  terms  of  the  autocorrelation  of  the  field  at  the  receiver  plane.  Goodman 
2 

has  shown,  in  the  case  where  the  real  and  imaginary  components  of  the  com- 
plex electric  field  obey  Gaussian  statistics  and  have  equal  variances  (which, 
for  example,  is  the  case  for  a diffuse  target),  that  the  speckle-intensity 
normalized  autocorrelation  function  [Tj(Ax,  Ay)]  (the  autocorrelation  co- 
efficient) and  normalized  field  autocorrelation  function  [r_(Ax,  Ay)]  are 

h 

related  by  the  following  expression: 

T_ (Ax,  Ay)  = <I(x,  y) ><I (x  + Ax,  y + Ay)>  + |T  (Ax,  Ay) | 2 
1 E (D-l) 

Taking  the  Fourier  transform  of  both  sides  of  this  expression  gives  the 
following  result: 


*J.  D.  Malick  et  al.,  Implications  of  Evolving  Laser  Technology,  Phase  II: 
Observables  and  Discriminants,  SAMS0-TR-71-220 , July  1971. 

2 

J.  W.  Goodman,  "Statistical  Properties  of  Laser  Speckle  Patterns,"  in 
Topics  In  Applied  Physics,  ed.,  J.  C.  Dainty,  Springer-Verlag : New  York, 

Vol.  9,  1975. 
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Figure  D-l.  The  Autocorrelation  Theorem.  This  theorem  (also  known  as  the  Wiener-Khintchine  Theorem) 
is  a general  result  of  Fourier  transform  theory  and  can  be  used  to  relate  the  optical 
properties  of  a target  to  the  observables  in  the  receiver  plane. 


TARGET  ELECTRIC  TARGET  BRIGHT-  AUTOCORRELATION  OF 

FIELD  DISTRIBUTION  NESS  DISTRIBUTION  TARGET  BRIGHTNESS 


sI(fx*  fy)  « «(fx.  fy)  + rx(Au,  Av) 


(D-2) 


where  S (f  , f ) is  the  power  spectral  density  of  the  speckle  pattern  and 
xx  y 

rj(Au,  Av)  is  the  autocorrelation  of  the  target  brightness  distribution, 
l(u,  v) . With  proper  normalization,  the  complete  relationship  is  as 
follows:* 


//l(  u. 


ST(f  , f ) = <I>  <6{f  , f ) + 
I x y | x y 


v) I (u  - XRf  )(v  - XRf  ) du  dv 
x y 


w 

JJ I(u,  v) 


du  dv 


(D-3) 

where  the  receiver  plane  spatial  frequencies  in  the  x and  y directions 

respectively  are  f and  f , the  mean  intensity  at  the  receiving-aperture 
x y 

plane  is  <I>,  the  perpendicular  distance  between  the  target  and  receiver 
plane  is  R,  and  the  wavelength  is  X.  The  second  term  in  the  braces  is  the 
normalized  autocorrelation  function  (the  autocorrelation  coefficient)  of 
the  target  brightness.  It  will  be  therefore  noticed  that  the  total  power, 
or  the  integral  of  the  entire  expression  over  frequency,  contains  equal 
contributions  from  the  DC  (impulse)  term  and  from  the  frequency-distributed 
component. 


With  this  relationship  between  the  speckle  PSD  and  target  brightness, 
it  is  straightforward  to  calculate  the  target  modulation  spectrum  for 
virtually  any  target  geometry  and  adaptive  illumination.  The  next  step  is 
to  relate  this  PSD  to  the  case  when  it  is  sampled  by  a finite  size  (rather 
than  a point)  receiver  and  how  this  is  scaled  into  the  temporal  frequency 
domain  to  account  for  target  rotation.  Only  the  final  relationships  are 
presented  here  and  the  reader  is  referred  to  Radley  et  al.*  for  the  de- 
tailed analysis. 


The  effects  of  aperture  averaging  on  speckle  are  shown  graphically 
in  Fig.  D-3  and  are  summarized  below.  The  power  spectral  density 

*J.  Radley  et  al. , COAT-Target  Interactive  Modulation  Effects.  GRC  Report 
571-02-TR,  August  1976. 


101 


SPATIAL  PSD,  6 (f 

9 n ' 


Figure  D-3.  Summary  of  Aperture  Averaging 
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S (f  , f ) of  the  speckle  intensity  incident  on  the  receiving  aperture  is 
ix  y 

related  to  the  target  brightness  distribution  I(u,  v)  and  its  mean  <I>  by 


SI(fx’  V = <X>  [6(fx*  V + ri(XRfx’  XRfy)]  (D"4) 

where  r_(ARf  , XRf  ) is  the  autocorrelation  coefficient  of  the  target 
i x y 

brightness,  X is  the  wavelength,  and  R is  the  range.  The  spatial  frequen- 
cies f and  f correspond  to  the  x and  y directions,  respectively,  in  the 
x y 

aperture  plane. 


The  power  spectral  density  S^C^,  fy)  of  the  speckle  intensity  that 
is  passed  by  the  aperture  is  related  to  S^(fx>  f^)  by  the  aperture  trans- 
fer function  H(f  , f ) 
x’  y 


vv  v - VV  VlH(V  v1 


(D-5) 


where  H(f  , f ) is  the  two  dimensional  spatial  Fourier  transform  of  the 
x y 

pupil  function  of  the  receiver  aperture.  The  power  spectral  density  (PSD) 
of  the  aperture  response  function  A(x,  y)  is  the  square  of  the  absolute 
magnitude  of  the  Fourier  transform  of  A(x,  y) . A(x,  y)  is  unity  within 
the  aperture,  zero  elsewhere. 

The  probability  density  function  p(I)  of  the  speckle  intensity  I at 
the  aperture  can  be  found  from  the  number  M of  correlation  cells  within  the 
aperture 


p(I)  = (MI/<I>)M  e-MI/<I>/ir(M) 


(D-6) 


where  F(M)  is  the  gamma  function  of  M,  and  M is  related  to  the  variance 
cf  the  aperture  integrated  speckle  PSD: 


M'1  ■ //Vf*’  V df*  df> 
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The  speckle  contrast  ratio  Cr>  defined  as  the  ratio  of  the  RMS  inten- 
sity variation  divided  by  the  mean  intensity,  is  related  to  M very  simply:  . 

C = M_1/2  (D-8) 

r 

When  there  is  relative  motion  between  target  and  aperture  (e.g.,  a 
rotating  target),  the  speckle  pattern  is  swept  across  the  receiving  aperture 
in  a direction  perpendicular  to  the  target's  effective  rotation  axis,  re- 
sulting in  a temporal-domain  PSD  SQ(f)  that  is  a function  of  the  temporal 
frequency  f,  and  is  given  by  a one-dimensional  integral  of  sq(^x»  f )• 

For  rotation  about  the  y axis,  the  expression  for  Sg(f)  is 


s0(f) 


(v  ) 

X 


-1 


/ 


S_(f  , f ) df 
0 x y y 


f — f /v 
x x 


(D-9) 


where  v^  = 2Rft  is  the  relative  velocity  in  the  x direction  at  the  aperture 
and  it  is  the  target's  effective  rotation  rate  perpendicular  to  the  line  of 
sight. 
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